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FUTURE USES OF LASER TRACKING 


E. M. Gaposchkin 
1. INTRODUCTION 

Laser technology has advanced to the extent that specific plans for its 
use can be made. Satellite tracking with lasers has become routine and the 
limits of accuracy are known, though only approximately realized. There 
are three different types of laser tracking units; the data are compatible and 
can be satisfactorily combined. 

Satellite tracking by laser is a very precise method of range measure- 
ment. The energy requirements are such that satellites must be equipped 
with corner-cube reflectors, which are passive devices and have an indefinite 
lifetime; there is no limit to the number of observations that can be made. 

Precise satellite tracking has evolved through many stages. Laser 
tracking is the most promising of the many ground-based techniques under 
development and in use today. The ultimate accuracy limit of laser tracking 
is set by the indeterminacy of the tropospheric refraction correction, a 
problem common to all ground-based observing systems. 

It seems, therefore, that we are at the final stage of ground-based 
tracking, and further advances will very likely come from satellite-based 
tracking systems. This discussion is limited to ground-based precision 
tracking and its viability. However, satellite -bas ed tracking will need to be 
related to the earth's surface. Precision tracking has been most commonly 
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used Ln geodetic studies. The first question is, thus, what impact will laser 
tracking have on determining the gravity field of the earth and the locations 
of observing stations. 

The present accuracy of laser ranges is 50 cm (see Section 5). Not all 
the systems in use have this capability, but even the primitive systems have 
an accuracy of 1.5 m. This can be compared with the 20-m accuracy of 
Baker -Nunn camera observations and the 15-m accuracy of the Goddard Range 
and Range Rate (GRARR) data. The ultimate accuracy of laser measurements 
will be 6 to 8 cm. To obtain this precision, we require very detailed mete- 
orological data, which are not yet available. 

There are many other uses for precision tracking of 50-cm accuracy. 

For example, ranging to the moon has been successful, and there is no doubt 
that uses other than those mentioned here will be found for laser satellite 

tracking. 


2 . SATELLITE GRAVITY-FIELD DETERMINATION BY CLASSICAL 

METHODS 

The methods discussed here are now classical (Kaula, 1966; Lundquist 
and Veis, 1966). The new data can be treated in the traditional framework, 
and alternative methods will be discussed in Section 7. 

2 , 1 The method used to determine the earth's gravity field is the inverse 
of the classical celestial mechanics problem, which is to determine the tra- 
jectory of a body, given a set of forces. -Observations are used to check the 
theory. The geodetic problem is to determine the forces, i. e. , the gravity 
field of the earth, given an orbital theory and observations of the trajectory. 
At the same time, one can and must determine the locations of the observing 
s rations . 
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2. 2 The representation of the geopotential in spherical harmonics, i. e. , 
in harmonic functions, is a spectral decomposition of the field. Specific 
features in the gravity field are not related to specific harmonic terms. 
Conversely, specific harmonics cannot be related to particular features. 

The situation is completely analogous to Fourier decomposition in one dimen- 
sion. In general, the harmonic functions that arise in the solution of the 
equations of motion are similarly not related to specific features in the 
motion, in this case a spectral decomposition in time. There are some pro- 
nounced spectral features in the motion, as we shall see, and hence to some 
extent the satellite acts as a filter, amplifying some spectral components. 

2. 3 We can determine the trajectory directly by numerical integration 
of the equations of motion. Alternatively, we can seek an analytical solution. 
In this case, we are forced to use approximate solutions, which appear to be 
satisfactory. This discussion is based on one such analytical theory. 


The zonal harmonics have the largest effects, and rather elaborate 
analysis is required. The equations of motion for a satellite under the 
influence of tesseral harmonics of the gravity field are approximated by 
expansions in harmonic functions. These expressions can be integrated 
directly to first order, which is satisfactory. 


2. 4 The equations of motion are solved by successive approximation. 
The expressions are for small changes, called perturbations, with respect 
to some reference, called a mean orbit. There are two classes of changes; 
1) nonperiodic (6 £ ), which increases indefinitely with time, and 2) periodic 


Periodic perturbations result from the spectral decomposition of the 
solution and are defined to have an average value of zero: 



(I) 
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The remaining motion is generally defined as the secular part. Terms 

? 3 

in t , t , etc. are not excluded in principle and apparently arise in practice. 

It is an open theoretical question if such terms can arise. They can also 

appear from very long-periodic terms that have not been included. In addition, 

some analytical theories have so-called mixed secular terms, e.g., t(sin at). 

For our purposes, we can ignore these questions and rely on equation (1) to 

define periodic perturbations, where t and t, are the first and last times of 

3. D 

the trajectory. Hence, if 6 £ is the total perturbation. 


o £ = 6 £ + 6 £ = 6 £'t + 6 £ 

s p s p 


then 


6 £ ' 
s 




a 


( 2 ) 


The 6£ g are integrated effects and become very large. Thus, the 

coefficients controlling 6£ , the zonal harmonics of even degree, can be 

s 

determined very precisely. These are integrated effects (in time) of specific 

harmonics (i. e. , averaged in space). For 6 £ , we can rely on integration 

P 

in time only for long-period terms. The latter fall into two categories: 1) 

zonal harmonics of odd degree, and 2) harmonics that are resonant with the 
satellite orbit. For the remaining terms in the harmonic decomposition of 
the trajectory, we cannot use integrated effects and must use the large spatial 
averaging of the harmonics, another kind of integrated effect. 

In summary, the gravity field is expressed in a spatial spectral decom- 
position, and the orbital theory in a temporal spectral decomposition. The 
satellite is sensitive to certain spatial components of the gravity field through 
certain temporal components arising in the orbital theory. 
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2. 5 The development of perturbations, now classical, is given in many 
papers (Kaula, 1966; Lundquist and Veis, 1966). The choice of dependent 
variables is arbitrary. The most common are to (perigee), £2 (argument of 
the node), I (inclination), e (eccentricity), M (mean anomaly), and a (semimajor 
axis). These elements can be combined into any other set. For descrip- 
tive purposes we use M and the shift in position dr = (dr • dr)^ = dr • dr. 


The expression of a periodic perturbation (say in M) due to harmonic 
is given in the form 


l 00 


6M 


lm = C lm 23 23 ' m ' Ps q> a> e ’ ^ - 2p) w + (l - 2p + q) M + m(Q - 0)] 

p= 0 q= - oo 


( 3 ) 


or 


6M. = C. 6m".' C[(i - 2p) w + (i - 2p + q)M + m(J2 - 0)1 

impq Im ImpqiJ 1 ' v /J 


( 4 ) 


where ^ is either sin or cos and 9 is the sidereal angle. 


A few well-known remarks follow: 

1. A(f , m, p, q, a, e, I) oc e ^ ^ ; hence the large st terms generally come 
with q = 0. 

2. The q summation, though formally from -coto +oo, only needs to go 
f r om -10 to +10. 

3. n = M > 9 = 1 » co, £2 ; therefore, the frequency 

f - (l - 2p)cu + (l - 2p + q)n+ m(£2 - 9) 

2 rr 

of any term is mainly controlled by 1VE and m9. 

4. Since A(£ , m, p, q, a, e, I) depends on elements that have virtually no 
change, A is constant for any particular satellite. The £ term contains all 
the temporal variations. 


5 


7 



2, 6 For any given value of order m, all perturbations of degree £ (even) 
will have the same frequency. Since 0 < i < m and 0 < p < 1 , arguments 
with f - 2p = r exist. This can be seen by examining the dominant term for 
each l, i.e., q = 0. The frequency becomes 

_ (l - 2p)(io + n) + m(fl - 9) 

4 _ 2tt 

With i = 12, p = 6, m = 12, we have the same frequency as with £ = 14, 
p = 7, m = 12; i = 16, p = 6, m = 12; etc. Similarly for £ (odd). In addition, 
we cannot have the same period for both £ odd and £ even. 

For example, the perturbations in M for DID (6701401) are given in 
expression (5) for only the principal terms with m = 1 , 2 and £ = 3, 4, 5, 6, 7, 8. 
For this satellite, a = 7614 km, e = 0. 0843, and I = 39.°455. 


6M = C 3 1 f - 7 - 1 sin (u + 8 - 0) + 0. 8 sin (to + 2M + 8 - 9) - 63. 3 sin (-gj + 8 - 0) + . . . ] 

+ Cj 2 ( -42 . 5 cos [to + 2 (£3 - 9)] + 1 0. 5 cos [go + 2M + 2(8 - 9)] - 13. 6 cos [-to + 2(8 - 0)] + . . . } 

+ C 41 [7. 0 cos (-M + 8 - 9) - 8. 2 cos (M + A - 9) + 5. 1 cos (-2w + 8 - 0) + . . . ] 

+ C 42 (-10. 3 sin [-M + 2(8 - 9)] + 14. 2 sin [M + 2(8 - 0)] + . . .} 

+ C 5J [-87. 4 sin (to + 8 - 9) + 6. 9 sin (to + ZU + 8 - 0) + 87. 9 sin (-w + 8 - 0) + . . .] 

+ C 52 (8. 6 cos jb + 2(8 - 0)] -1.4 cos [w + 2Mf 2(8 - 0)] f 43. 9 cos [-oj + 2(8 - 0)] + . . . } 

+ j [ 5. 1 cos (-M + 8 - 9) - 6. 0 cos (M + 8 - 0) - 16. 2 cos (-2w + 8 - 0) + . . . ] 

+ C 62 {5. 4 sin [ - M + 2(8 - 0)] - 7. 4 sin [M + 2(8 - 0)] + . . . } 

+ C_, j [ 33. 1 sin(co + 8-0) + 0. 0 sin (u + 2M + 8 - 9) + 1.4 sin (-to + 8-0)+...] 

+ C^ 2 {40. 0 cos [to + 2(8-0)] - 5.5cos[to+2M+2(8-9)] - 40. 3 cos[-oo+ 2(8- 0)] + ...} 
+ Cgj[-6.8 cos (-M+8-0J+ 7.9 cos(M + 8-0) + 19.1 cos (-2o + 8 - 0) + . . . ] 

+ C g2 {4. 1 sin [-M + 2(8 - 0)] -5.7 sin [M + 2(8 - 0)] + . . . } 

+... . ( 5 ) 
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We can rearrange this expression in terms of the same frequency (the period 
of each term, in days, is given in parentheses): 

6M = (-7. 1 C 31 - 87. 4 C gl + 33. 1 C yl + . . . ) sin (go + £2 - 0) (-1. 001 ) 

+ (0. 8 C 31 +6.9 C 51 + 0. 0 C ?1 + . . . ) sin (go + 2M + £2 - 9) (0. 040) 

+ (-63.3 + 87. 9 + 1.4 C y ^ + ..) sin (-go + £2 - 9) (-0. 971) 

+ (7. 0 + 5.1 - 6.8 Cgj + ...) cos (-M + £2-9) (-0. 071) 

+ (-8. 2C 41 - 6. 0 C 6I + 7. 9 C + . . . ) cos (M + £2 - 9) (9. 083) 

+ (5.1 - 16.2 + 19. 1 Cg^ + ...) cos (-2 go + £2 - 9) (-0. 958) 

+ (-42. 5 C 32 + 8. 6 C 52 + 49. 0 C ?2 + . . . ) cos [go + 2(£2 - 9)] (-0. 49?) 

+ (1 0. 5 C 32 -1.4 C 52 - 5. 5 C ?2 + . . . ) cos [w + 2M + 2 (£2 - 9)] (0. 041) 

+ (-13. 6 C 32 + 43. 9 C 52 - 40. 3 C ?2 + ...)cos f-u+ 2(£2 - 9)] (-0. 327) 

+ (-10. 3 C 42 + 5. 4 C 62 + 4. 1 C g2 + . . . ) sin [-M + 2(£2 - 9)] (-0. 066) 

+ (14. 2 C 42 - 7. 4 C 62 - 5. 7 C g2 + . . . ) sin [M + 2(£2 - 9)] (0. 091 ) 

+ . . . . ( 6 ) 

As a second example, we give in expression (7) the perturbations in M 

✓ 

for the satellite PEOLE, to be launched by Centre National d'Etudes Spatiales 
(CNES), for which a = 7173 km, e = 0.009, and I = 15.°0. In this case, we 
give the terms for m = 1 , 2 and l = 3, 4, 5, 6. 
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OM - C ^ j [-1493 sin (u +32 - 6) + 160 sin (to + 2M + 12 - 0) - 164 sin (-<o + 12 - 0) + 

+ C 32 j-341 cos [<o + 2 (i2 - 0)] + 77 cos [ u + 2M 1 2 (12 - 0)] + . . . | 

+ [126 cos (-M + 12 - 0) - 146 cos (M + 12 - 0) + . . . ] 

+ C 42 | 1 35 sin [2io + 3M + 2(12 - 0)] - 35 sin [-M + 2(12 - 0)] + 47 sin [M + 2(<> - 0)] 

+ 50 sin [ 2io + M + 2(12 - 0)] + . . . | 

+ C,-j [1600 sin (to +12 - 0) - 114 sin (go + 2M + 12 - 0) + 579 s in (-go + 12 - 0) + . . .] 

+ C' 52 1 761 cos [w + 2(12 - 0)] - 115 cos [go + 2M + 2(12 - 0)] + . . . | 

t C 2 , [-165 cos (-M +12 - 0) + 190 cos (M +12 - G) + . . . ] 

t C 62 | -7 6 sin [2 + 3M + 2(1 2 - 0)] + 73 sin [-M + 2(12 - 0)] - 99 sin [M + 2(12 - 0)] 

- 69 sin [2co + M + 2(12 - 0 )] + ... | 

(7) 

We can rearrange this expression in terms of the same frequency (the period 


of each term, in days, is given in parentheses): 

6M = (-1493 C 3 + 1600 C 5 + . . . ) sin (cj + £2 - 0) (-1.013) 

+ (160 C 31 - 1 14 C 51 + . . . ) sin (w + 2M + £2 - 0) (0. 0 36) 

+ (-1 64 + 579 + . . . ) sin (-co + Q - 0) (-0.949) 

+ (126 - 1 65 + ...) cos (-M + 12 - 0) (-0.065) 

+ (-146 + 190 + ...) cos (M + i~2 - 8) (0.07 5) 

+ (-341C 32 + 76lC 52 + ...) cos [oj + 2(12 - 0)] (-0. 498) 

+ (77 C 32 - 115 C 52 + . . . ) cos [ w + 2M + 2(12 - 0)] (0.038) 

+ (135 C^ 2 - 76 C^ 2 + ...) s in [ 2u> + 3M + 2(12 - 0)] (0.024) 


10 


8 



(- 0 . 061 ) 


+ (-35 C 42 + 73 C 62 + . . . ) sin [ -M + 2 (ft - 0)] 

+ (47 C 42 - 99 C 62 + . . . ) sin [M + 2 (ft - 0)] (0.082) 

+ (50 C 42 - 69 C 62 + . . . ) sin [2co + M + 2 (ft - 0)] (0. 081) 

+ . . . . ( 8 ) 


Even if we assume the satellite to be a perfect filter, uncontaminated by 
other model errors, and the tracking data and analysis process to be perfect, 
we see that with one satellite we can only determine spectral components that 
are linear combinations of the gravity field (C^ ) and functions of orbital 

elements [A(f , m, p, q, a, e, I)] . From each satellite we obtain one or two 
linear combinations of harmonics for l odd and for i even. By using addi- 
tional data we can only refine the numerical value of these linear combina- 
tions. The coefficients of the relations will depend on the orbital elements 
so that other linear combinations can be determined only from additional 
distinct orbits. Generally, this is achieved by selecting satellites with 
different inclination, but independent linear relations can also be obtained 
with changes in eccentricity e or semimajor axis a. 

As the degree increases, the perturbations become negligible, and so the 
linear relation does not involve an infinite number of parameters. The num- 
ber of parameters required is determined from this study. 


From equation ( 6 ) (for DID) we see that a linear combination of C,,, 

3 i 

C 51 , Cyj , ... can be determined from the - 1.001 -day-period term and 
another of equal size from the -0.971 -day term. The third term (and there 
are many similar smaller terms) is a factor of 10 smaller and will not con- 
tribute significantly as an observation equation. The linear combination of 


c 32 , c 52 , c ?2 , . . . 

-0. 327-day period. 


has only one significant spectral component for the 


9 


11 



F rom equation (8) (for PEOLE) the linear combination will be determined 
from the -1.013-day period. There are two other components at periods of 
0.036 and -0.949 day. The effects are a factor of 2 to 10 smaller. A second 
linear relation of lower weight can be determined. These additional terms 
arise in satellites with nonzero eccentricity; e.g., for DID the terms of 
q + 0 are significant, as illustrated by equation (6). 

Therefore, one or two finite linear relations are determined for £ (odd) 
and one or two for £ (even). In addition, weaker relations can be established. 
Each satellite can contribute to the unique determination of 1 or 2 odd and 
1 or 2 even degree harmonic coefficients in each order; i. e. , if there are 
24 sets of unique gravity-field coefficients affecting the orbits for a given 
order, then 12 distinct satellites would be sufficient to determine them. 

For a nonhomogeneous set of satellites, i. e. , where they are not all 
equally sensitive to the gravity field, subsets of coefficients are determined. 
In'the case where insufficient satellites are available, the linear relations 
are generally solved by constraining the higher degree and order coefficients 

to zero. 

The linear relations are not determined with equal accuracy; for exam- 
ple, the resonant harmonics have a very large effect and the spectral compo- 
nent is strongly determined. However, the resonant period is commensurate 
with the arc length, which will cover only a small number of cycles. This 
makes separating nearly commensurate periods difficult. For spectral 
components bunched between P = 0. 02 and P = 0. 04 day, i. e. , between 50 
and 6 0 min, the effects are small and the spectral decomposition is also 
difficult. 

3. GRAVITY-FIELD DETERMINATION BY TERRESTRIAL MEASUREMENTS 

3 1 The gravity field of the earth is measured with gravimeters. These 
measurements have been made on part of the surface only — 25% of the 1° X 1° 
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elements of the earth's surface. By combining these into 5° X 5° elements, 

56% of the surface is covered. These measurements can be used to deter- 
mine the spherical harmonic expansion representing the gravity field. 

3 ■ 2 Rigorous analysis of surface -gravity measurements can lead to the 
determination of these spherical harmonic coefficients. This requires 1) 
knowledge of gravity everywhere on the surface (e. g. , with 5° X 5° averages), 

2) reduction of gravity values to a common reference, and 3) knowledge of 
the correlation function (now generally approximated) (Heiskanen and Moritz, 
1967). 

3 . 3 The above three conditions are not met today, so the usefulness of 
surface-gravity measurements is reduced. For example, one cannot deter- 
mine spherical harmonics, from the data available, without some assumptions 
such as model anomalies. Used in combination with other techniques, surface- 
gravity measurements provide additional observation equations and comparisons. 

3 . 4 Local gravity cannot give much information about very long-wave- 
length properties of the gravity field. The lower degrees £ = 2, 3, 4, 5 are 
poorly determined; they are, however, particularly well determined by satel- 
lite techniques. Apart from this consideration, there is no preferential set 
of harmonics determined by surface gravity, unless there is a "topographic" 
selection. The satellite effects for degrees £ = 10, 11,. ..,25 are dominated 
by resonances, and there is virtually no information for C/S (f, 5), C/S (£, 6), 
and C/S (i, 7) for i = 10, 1 1 , . . ., 25. 

3. 5 The combination and comparison of surface-gravity anomalies 6g 
and geopotential anomalies 6ifi involve the gradient of tujj: 

6g = y(i + 1) 6 i|j . ( 9 ) 

Hence, the variance of 6 lJj, determined from 6g, increases more slowly than 
the variance of 6g by the factor 1 / (£ + 1). The effects of otjj from satellite 
perturbations become very small as £ increases. Therefore, the lower 
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harmonics will be determined from satellite studies, and the higher harmonics 
from surface gravity (Section 3.4). 

3. 6 Rapp (1968) has pointed out that in surface-gravity analysis the 
amplitude of the spherical harmonic coefficient decreases with increasing i, 
the variance increases with increasing l, and they cross at approximately 
i = 18. If one adopts the position that the coefficients cease to be meaningful 
when they are smaller than their variances, then f = 18 seems to be the limit 
of usefulness of surface-gravity measurements. 

3 . 7 The rule of thumb for decrease of coefficient size is 


C 


i 


m 



( 10 ) 


Hence , C ,, =3X10 ® which is equivalent to 980 X (18 + 1) X 3 X 10 ^=0.56 

1 8, m 

mgal (from equation (9)). If we consider the result for the variance of a 

5° X 5° square. 


2 1000 .2 

cr = mgal 

g n 


where n is the number of anomalies given in the 5° X 5° square, then for 

n = 25, we have 


cr = 's/I 000/ 25 = 6. 3 mgal 

g 

Even if we take the optimistic value of 6. 3/4 = 1.6 mgal, it hardly approaches 
the 0.56 mgal required. 
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4. 1969 SMITHSONIAN STANDARD EARTH (SEII) 


4. 1 The SEII (Gaposchkin and Lambeck, 1970) was determined primarily 
from optical tracking data, combined with significant but not dominating 
amounts of laser data. The number of optical data used in the form, of simul- 
taneous observations was comparable to that used in a dynamical solution. 

The former give information about relative station positions, and the analysis 
is relatively free of assumptions. 

4. 2 The determination of station coordinates is significantly improved 
by the use of both geometrical and dynamical methods for two reasons. First, 
a completely different analysis for simultaneous observations was not con- 
taminated by assumptions made in the dynamical solution. For example, 
some orbits were chosen because of large gravity-field effects. These effects 
degraded the solution for coordinates during early iterations and the simul- 
taneous observations canceled out these biases. Second, the validation of the 
results is critical and was made possible by the two methods. 

4. 3 A further geometrical data set was included — the relative longitudes 
and distances to the axis of rotation of the JPL DSN stations — which strength- 
ened the solution and provided an additional check on the results. 

4. 4 Twenty-one satellites were used in the dynamical analysis. Only 
six of these had any laser data. In addition, there were only five laser 
stations. The 21 satellites ranged from 28° to 95° inclination. The variation 
in orbital elements was sufficient for the separation of the harmonics. The 
complexity of the field determined was fixed by the accuracy of the data (15 
to 20 m). The gravity field could be completely determined through £ = 10, 
m = 10 from satellite data. For the 11th- and 12th-degree terms only 
C/S (11,7), C/S (12,6), and C/S (12,9) we re not determined. Highe r o rder 
terms selected were C/S (1,1), 1 3 < £ < l6;C/S(i,2), 1 3 < I < 15; C/S 
(14, 3); C/S (i, 12), 13 <i < 19; C/S (i, 13), 13 <f < 21 ; and C/S (i, 14), 

14 :S i < 22. The 12th-, 1 3th-, and 14th-order terms are, of course, the 
resonant harmonics. 
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4, 5 The remaining harmonics through l - 16, m = l6 are determined 
in a combination with the available surface -gravity data: 935 300 X 300 n mi 
squares out of the 1654 possible squares. Of these, only 136 had more than 
20 T X 1° anomalies. The 1 6th degree and order were chosen to be the point 
where the data did not significantly contribute to an extended field; i. e. , the 
residuals were not reduced by adding unknowns. In addition, some surface- 
gravity anomalies were rejected on a 3cr criterion. Generally, these rejected 
anomalies were in regions of significant tectonic activity, such as the Puerto 
Rico trench. Hence, the rejected data are probably good, but not character- 
istic of the long- wavelength features of the gravity field that we are deter- 
mining. 

4. 6 Figure 1 shows the degree variance and the standard error of the 
coefficients for SEII. The coefficients for i > 16 are determined by the 
resonant effects. The standard error of the solution for the harmonics will 
apparently cross the degree variance line between 1=18 and I = 20, con- 
firming Section 3.6. Since the harmonics from i = 1 2 to f =16 were mostly 
determined by surface gravity, this result is not independent of Section 3.6. 

5. LASER TRACKING 

5 . .1 The precise measurement of range gives one component of the satel- 
lite position. For zenith passes this will be a radial component dr that is 
less affected by the earth's gravity field. Observing at lower elevation angles 
measures part of the along-tracking component du, which is generally the 
largest effect, and the across-track component dw, which is the same size as 
dr. For the following discussion, we assume that we measure a general 
component of the position and are not restricted to the radial component. If 
we observe at an elevation of 30°, we get \l 3/2 =0.87 of the along-track 
component, 

5 . 2 Laser systems differ in three principal ways: 

A. Accuracy. The first systems were limited to an accuracy of 1 . 5 m, 
and some are still at that level. Newer systems have achieved an accuracy 
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of 50 cm (Lehr, Pearlman, and Scotty 1970; Pearlman, Lehr, Mendes, and 
Wolf, 1970). The limitation on accuracy arises from two sources. The first 
is the resolution of the equipment itself and is discussed in Lehr et al . (1970) 
and Pearlman et al . (1970). The second, and ultimately limiting, source of 
error is the determination of refraction for the atmosphere. The current 
accuracy (Lehr et al . , 1970; Pearlman et al ., 1970; Lehr, Pearlman, Salisbury, 
and Butler, 1969) of the refraction correction is sufficient. For future observa- 
tions where system accuracy is improved to between 1 and 5 cm, the atmos- 
pheric correction will have to be treated with greater precision, perhaps 
along the lines described in Hopfield (1970). 

B. Repetition rate. The number of observations per pass run from 5 
to 200. The passes with 200 observations have been analyzed statistically, 
and the observation errors (residuals) have a white noise spectrum. 

C. Acquisition method. The SAO and Goddard Space Flight Center 
(GSFC) systems have automatic - pointing capability and can observe 24 hours 
a day. Other systems, e. g. , those operated by CNES and the National 
Technical University (NTU), Athens, acquire the satellite visually and can 
observe only during twilight. 

Some of these differences are related. For example, the number of 
observations is related to the cooling system, as is the power output. The 
latter determines the efficacy of daylight observations and accuracy. 

5 . 3 The minimum elevation for successful observations could be as low 
as 15° or as high as 30°. The ability to make the atmospheric correction is 
the controlling factor. 

5 . 4 In order to use laser tracking data to determine the geopotential, 
we must sample the gravity field everywhere. We can find out how many 
stations are necessary for complete coverage. Table 1 lists for a variety of 
heights (h) the geocentric angle subtended (9), the fraction of the earth's 
surface sampled by one station (£,), and the number of stations necessary for 
an almost complete sampling of the whole earth (N ). Each quantity is given 

S L3. 
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for zenith distances of 60° and 75°, corresponding to elevations of 30° and 
15°, respectively. Also computed is the subsatellite height at minimum 
elevation (h'). The satellite must reach this height to be observed every- 
where at minimum elevation. Therefore, h' would correspond to a minimum 

apogee height. 

5 5 Table 1 is computed for a polar satellite, and Figure 2 represents 
the computation graphically. For a satellite of arbitrary inclination the 
amount of the earth covered will be less. These figures can be appropriately 
modified by multiplying [ t, sin (I + 0)] by [N^^/sin (I + G)]. 

5.6 It is quite apparent that to obtain reasonable coverage with a modest 
number of stations (say, <50) for satellites with low apogee heights (say, 

< 1 Mm), we must make observations at the minimum elevation of 15°. 

5 . 7 This analysis holds for any ground-based tracking system. 

6. FUTURE GRAVITY -FIELD DETERMINATIONS WITH 

LASER SATELLITES 

6. 1 In Section 2.4 the decomposition of a satellite perturbation into a 
sum of products was described. Each term in the sum is of the form 

6 £ . =S„ 6 £* + pm + - 6)] , 

Impq lm lm v L J 

where S „ is the gravity-field coefficient; 6 £ is a function of a, e, and 

lm 6 1 impq 

I and is virtually constant for any satellite; and/Tfaco + Pm + y(S2 - 0)] is a 
trigonometric sine or cosine function of the angular variables, a,p,y being 
integers. Since -1 < 1 and the amplitude of S^ m can be roughly esti- 

mated as 


S, oc 10 5 /| 2 
lm 
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we can estimate which terms 6£„ are significant by calculating 6£. 

|mpq t, i t. impq 

Any C/S (i, m) may be much smaller than the estimate, and the perturbation 

will not, in fact, be necessary. Conversely, there could be very large 

- 5 2 

C/S (i,m), which would be ignored by use of the rule (10 /i ). These two 
possibilities must be reinvestigated after the coefficients are determined 
by examining the numerical values or by looking for unmodeled effects in 
the orbital residuals. 

The term 6£„ is a generic element, and it is more meaningful to 

impq b b 

collect the perturbations of the elements into the perturbation of the position 

vector with components du, dw, and dr. The total change is computed as 

(du^ + dw^ + drV ^ = dp. Values of dp X 10 ^/l^‘ ^ >0.01 X 10 ^ are kept 

1.7 

as significant. The value of I is used as a conservative estimate of the 
decrease in the size of C/S (i,m). We therefore tend to keep more coefficients 
than necessary. For future reference, we give in Table 2 the estimated size 
of C/S (i , m) . 

In Tables 3 to 7 and 9 to 1 2 the 5 £ . are collected into dp and summed 

impq r 

over p and q for any i,m. The units are given so that the product of the 

6 

numerical and gravity-field coefficients multiplied by 10 yields the pertur- 
bation in meters. For example, in Table 3 (PEOLE satellite) for 1 = 3, 
m = 1 , dp = 271, and C ^ ^ =1.9X10 ^ , the total perturbation due to C } is 
1.9 X 271 =515 m. From one aspect, the coefficients are sensitivity 
coefficients . 

6. 2 Table 3 gives the complete list of sensitivity coefficients for the 
proposed PEOLE satellite through i = 20, rn = 20. We note two characteristics 
of PEOLE: 1) for any degree i the size of the perturbation decreases with 

the order, and 2) there are no effects for m 2: 7 for any degree, but all m > 6 
are important through & = 20. 

Table 4 lists sensitivity coefficients for satellite 6508101 (OGO 2) for 
11 < i < 20. The main differences between OGO and PEOLE are shown as 
follows : 
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PEOLE 


OGO 


I 

a 

perigee 

apogee 


15 ° 

7. 173 Mm 
730 km 
860 km 


7.343 Mm 
424 km 
1 504 km 


Note that for OGO all but a few coefficients are significant, although it has 
a larger semimajor axis. This Ls essentially due to the high inclination, 
because OGO samples the field everywhere. However, the size of the lower 
order terms is smaller. Therefore, PEOLE will be stronger in determining 
coefficients m < 6. Also, OGO has resonant harmonics of the 14th order. 
OGO is not a laser satellite and is given here for illustrative purposes only. 


Higher inclination satellites sample more gravity-field coefficients, 
and lower inclinations help to separate the coefficients. 


6 . 3 Table 5 gives the sensitivity coefficients for satellite 6701401 (DID), 
a laser satellite of moderate height at I = 39°. Note that 1) it does not sample 
the complete field even to l - 10, m = 10; 2) most of the coefficients to i =15, 
m = 15 are sampled; and 3) between f = 17 and i = 20 a little more than half 
the coefficients are sampled, and even for i = 20, m = 5 the effect is 
0.025 X 18 = 0.45 m. 


Table 6 gives the coefficients for satellite 6503201 (BE-C). We can 
determine the complete field with laser data only to perhaps i = 15, m = 15. 

6 . 4 Table 7 gives the sensitivity coefficients for the seven laser satel- 
lites for l = 20. As we previously noted, not all the coefficients are sampled: 
m = 7, 9, 16, 17, 18, 19, 20 are not determined at all; m = 6, 8, 10 are very 
weakly determined. Apart from the resonant harmonics for l > 20, a 
reasonable field to I = 20, m = 20 would allow orbit calculation to an accuracy 
of the data. 
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6. 5 The seven satellites can determine 7x2 = 14 degrees of gravity- 

field coefficients. It is clear that the satellites at different inclinations are 

distinct (see Section 2. 5). We will show that the three satellites at 40° are 

also distinct. This can be seen by examining the ratio of the perturbation 

for each order. The analogy of spectrum analysis shows that each order 

provides a linear equation relating the coefficients of all degrees. If two 

satellites are identical, these two equations will be the same; i. e. , the 

coefficients of C, will be in the same ratio. We then investigate the ratio 
* irr i 

of 6p (DlD)/6p^ (BE-C) in Table 8. We see for the odd d egrees that 

x j m x y in 

the ratios are significantly different, as they are for the even degrees. 

Hence, the seven satellites do indeed constitute seven unique objects. 

6. 6 This selection of satellites will allow us to determine the field 
approximately to f = 18, m = 18. However, there are significant effects for 
higher order harmonics that are at best weakly determined. We would suggest 
one or more additional satellites to aid this analysis. A polar satellite would 
be the best choice since it samples the field everywhere. 

Tables 9 and 10 give the effects for the 18th-degree harmonics for an 
inclination of 85° and a variety of apogee and perigee heights. It is clear that 
the 450- by 500-km orbit has the largest sensitivity to the geopotential. 

Perigee was arbitrarily chosen at 450 km. The air drag for a geodetic 
satellite at that height can be adequately treated theoretically. For example, 
the change in semimajor axis for satellite 6701401 is less that 1 m per revo- 
lution, and an analytical theory can be used to compute that perturbation. 

The general "secular or logarithmic" terms are satisfactorily computed with 
empirical polynomials. 

For a complete sampling of the gravity field directly by use of geodetic 

satellites such as 6701401 to 18th degree and order, the apogee height must 

be 550 km. From Figure 2 , if d> = 7 5°, we would need about 120 stations, 

b T max 

clearly an unreasonable number. 
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6 , 7 Determination of spherical harmonics to i = 18, m = 18 describes 
the gravity field to a linear resolution of 

If X 6378 = 2220 km , 

A O 

which is equivalent to a central angle 0 = 19.° 5. Referring to Table 1, if 

4> = 75° , we need roughly 30 stations and satellites with apogee heights 

m ci x 

(corresponding to h') of at least 1. 3 Mm. 

7. ALTERNATIVE METHODS OF GRAVITY -FIELD DETERMINATION 

BY SATELLITE TECHNIQUES 

7. 1 Satellite altimetry to determine the sea surface to an accuracy of 

1 m and satellite-to-satellite tracking promise to provide significant improve- 
ments in our knowledge of the geopotential. Neither method yields informa- 
tion about site coordinates. Both methods have other uses, such as routine 
orbit maintenance. The satellite-to-satellite tracking used to improve the 
gravity field will need low satellites for the same reason that classical 
methods require lower satellites. The principal advantages of these methods 
are complete coverage and reduction of some refraction errors. 

7 . 2 These two methods of gravity determination will apparently become 
useful operational systems within the next 5 years. They will never improve 
knowledge of the long- wavelength phenomenon that is determined from long 
integration effects in the classical method. Therefore, the three methods 
should be combined. 

7 ■ 3 These two new methods require precise tracking from the ground, 
which, presumably, will be provided by laser and very long baseline inter- 
ferometry (VLBI) tracking. It is apparent that laser tracking will become 
a vital component in these future methods. 
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8. FURTHER USES FOR LASER TRACKING 


8. 1 Any program that uses precise positional information on a global 
scale will ultimately use laser tracking in some form. The GEOLE project 
(Husson, 1969) is one such program. Laser ranging, perhaps to the moon 
or to synchronous satellites, can measure the polar motion and rotation of the 
earth. Of course, laser tracking is ideal for routine orbit maintenance. The 
satellite system is completely passive. It provides day and night operation 
with reasonably simple equipment. 

8. 2 We now turn to the discussion of high satellites, which are very 
stable orbiting objects with minimal perturbations from the gravity field, 
air drag, and radiation pressure. Such satellites of intermediate height can 
be used for determining polar motion, variations in the rotation of the earth, 
station location, and ultimately the temporal motion of observing stations, 
which will measure crustal motion and earth tides. 

Our present knowledge of the gravity field limits the usefulness of such 
a satellite. Tables 11 and 12 list the sensitivity coefficients for two satellites 
of this type. As expected, the low-order terms are the largest. There are- 
modest and manageable resonances. Table 13 lists the formal uncertainties 
for the first four degrees resulting from the SEII solution. If we conserva- 
tively take 0.01 X 10 ^ as the uncertainty in any C„ /S. , then for SAOl we 

would make an error of 1.7 m for or 1.4m for ; for SAG2 this error 

reduces to 0.8 m for C It is apparent that improvements in the gravity 
field of the earth in the next year or two will easily reduce this uncertainty 
by a factor of 10, which allows an accurate enough determination of the 
gravitational perturbations to calculate the position of the satellite to 20 cm. 

The principal problem with a high satellite will be the radiation pressure. 
Even with extreme mass-to-area ratios, the long-term effect cannot be 
eliminated. However, with suitable satellite design this perturbation can be 
computed with sufficient accuracy. For example, for a satellite with 


21 


23 



M/A = 3300 kg/m^, I = 80°, e = 0.05, perigee = 2500 km, and a = 9.345 Mm, 
the short -period effects will be 6.5 cm in one revolution. The maximum 
perturbation, which is really a very long-period effect, will be 40 cm after 
one revolution. The short-period perturbation in the radius vector in one 
revolution will be less than 1 cm. This satellite would be ideal for studying 
the vertical component of the earth tide, which is 50 cm. 

Laser tracking of synchronous satellites would have the advantage of 
reducing the effects of the uncertainties in the earth's gravity field. The 
radiation-pressure effects are greater since they are proportional to a. 
However, the geometry of ranging is significantly degraded. For geometrical 
programs the heights of SAOl and SA02 are preferable. 

9. CONCLUSIONS 

A. Even with 10-cm laser observations the gravity field cannot be 
determined complete to £ = 20, m = 20. 

B. By use of optical data on lower satellites, in conjunction with laser 
data, the gravity field can be determined accurately enough to utilize fully 
IQ-crn laser data. With higher satellites the gravity field would be accurate 
enough to compute 15 -cm orbits. 

C. In addition to the resonant harmonics, many higher harmonics will 
have to be determined. This can be done by combining surface-gravity 
data, optical data from lower satellites, and laser data. 

D. Kinematic properties such as earth tides, polar motion, and crustal 
motion can be studied by means of these accurate laser data, relatively 
unaffected by lack of knowledge of the gravity field. 

E. Without large quantities of laser data the goals of the National 
Geodetic Satellite Program (NGSP) for a i = 15, m = 15 gravity field cannot 
be reached without the use of surface-gravity data. 
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F. Currently available surface -g ravity data (classified or not) used in 
combination with satellite data will not permit determination complete to 
£ = 20, m = 20. 


G. The NGSP goal of 10 m for station-coordinate accuracy has been 
reached for some stations. The ultimate accuracy of laser sites can be 20 
to 30 cm. For this accuracy to be achieved, an amount of data comparable 
to that used in the SEII must be acquired. 

H. To attain the i = 20, m = 20 gravity field, we must determine 420 

coefficients, requiring a computer memory of 90,000 words. In addition, 

to determine station coordinates and higher harmonics, say 30 stations and 

2 

60 additional harmonics, we require (420 + 30 X 3 + 60) /2 = 163,000 words 
of computer memory. 

I. If we use satellites with apogee heights of 1 . 3 Mm, we require 30 
stations to obtain complete global coverage of the gravity field, assuming no 
observations below 15° elevation angle. 

J. The station selection should be such as to distribute the area covered 
globally; the Southern Hemisphere in particular has been inadequately covered 
in the past. In addition, as many as possible of the original Baker-Nunn 
sites should be retained to allow better initial starting values and a tie to an 
inertial reference system. 

K. Assuming that the seven satellites described in Table 7 are distinct 
for purposes of gravity-field analysis, we can separate 14 harmonics. By 
use of optical observations on other satellites, the lower harmonics are 
separated. An advisable approach is to add three laser satellites at distinct 
inclinations, say at 30°, 50°, and 90°. These could be retrograde or pregrade. 
Retrograde satellites provide additional information on the odd zonal harmonics 
and a consistency check on timekeeping, both clock time and UT1. 
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L. For orbital maintenance, such as GEOLE, ATS-F, or ATS-G, fewer 
laser stations would be needed. Indeed, it seems that the geodetic program 
is making the greatest demands on laser satellite tracking. 
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Figure 1. Degree variances for 6 < i < 22 for the combination solution. 

Kaula's rule of thumb is indicated by the dashed line. For 2<i<6 
the degree variances are in complete agreement with this rule. 
The lower curve gives the degree variances corresponding to the 


precision estimates of the harmonics, i. e. , 1 /n' ^ + cr^ 
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Table 1. Dependence of station visibility on satellite height. 
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Table 2. Estimated size of C„ /S. = 10 ^/f 2 or 10 

im fm 

in units of 10”^. 


i 

10/f 2 

,0/i 1 - 7 

2 

2. 50 

3. 08 

3 

1.11 

1 . 54 

4 

0. 625 

0. 947 

5 

0. 400 

0. 648 

6 

0. 278 

0.476 

7 

0. 204 

0.366 

8 

0. 156 

0. 291 

9 

0. 123 

0. 239 

10 

0. 100 

0. 199 

1 1 
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0.169 

12 

0.069 

0. 146 
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0.059 

0. 128 

14 

0.051 

0.113 
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0. 044 

0. 100 
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0. 039 
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Table 3. Sensitivity coefficients for satellite PEOLE. 
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Table 4. Sensitivity coefficients for satellite 6508101 (OGO 2). 
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Table 5. Sensitivity coefficients for satellite 6701401 (DID). 
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Table 6. Sensitivity coefficients for satellite 6503201 (BE-C). 
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Table 7. Sensitivity coefficients for seven laser satellites for 20th-degree harmonics. 
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Table 9- Sensitivity coefficients for three hypothetical 
satellites of 85° inclination for 18th degree. 
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Table 10. Sensitivity coefficients for six hypothetical satellites of 85° inclination 
fori 8th degree. 
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indicates a very long-period term. 



Table 11. Sensitivity coefficients for satellite SAOl. 
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Table 12. Sensitivity coefficients for satellite SAG2. 
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Table 13. Formal statistics for SEII gravity field. 
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Abstract 


AFCRL has been actively engaged in the instrumentation 
of a system for laser-satellite geodesy. The system is a dual, 
laser concept to obtain range using a Q switched laser and to 
obtain angular information by photographic reflected high energy 
normal mode laser pulses against stellar fields. The Q switched 
laser system is capable of obtaining up to ten range measure- 
ments in a single pumping period which increases the confidence 
in returns, especially in the presence of high noise. Satellite 
range measurements have been made using only satellite position 
predictions. 
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AFCRL has been engaged in the development of a geodetic laser 
system in which the approach has expanded into a dual laser concept. 

One laser, operating in the Q switched mode, is used for satellite 
ranging. Another laser is operated in the normal mode for satellite 
illumination - sufficient for photographing the reflections in star fields. 
One advantage of such an approach is that fewer stations need to make 
observations for geodetic reductions. This means that for a given 
number of stations the time required for geodetic positioning is reduced, 
mainly because of weather restrictions, i. e. a higher percentage of 
stations can be "weathered in" without seriously affecting the time scale. 

Other aspects of the system differ from "normal" in other ways also. 
The ranging laser is capable of obtaining up to ten pulses in a single 
pumping period allowing ten range measurements to be made in 300 
microseconds. The Haute Provence dual laser system operates 
alternately - a range measurement is made and 4 seconds later a 
normal mode laser is fired for a photograph; 4 seconds later a range 
measurement, etc. - so as to avoid return echo ambiguity. It is felt 
that the AFCRL system will avoid ambiguity by virtue of the ten range 
measurements made by the multipulse laser. Therefore, the range and 
angle measurements can be made nearly simultaneously. 

The range laser is Q switched with a Pockels cell ten times per 
pumping period, normally at a repetition rate of once per 30 micro- 
seconds (the rate can be varied). The output is about one joule per 
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pulse with half power points of 30 ns. The output is presently collimated 
to two milliradians. When the cooling system is completed the laser can 
be pulsed once every three seconds. 

Start pulses are generated when the laser is fired. The pulses are 
picked off the rear of the ruby, sent to a photodiode through fiber optics. 
The return is collected by a 22 cm telescope and then directed to an EMI 
9558B photomultiplier after passing through a 10 A filter. The PMT high 
voltage is normally off and is gated on for five ms during the time the 
return is expected. A mechanical shutter is also used if the normal mode 
laser is operated immediately before or after the range laser but before 
the return is expected. It has been found that, even with the high voltage 
off during the time of firing a high energy laser, the PMT is extremely 
noisy when the high voltage is gated on up to 20 ms later. 

The elapsed transit times are measured by a counter system with 
one ns resolution. Since there are 10 expected returns at specific 
intervals the counter system must be range-gated 10 times. Range- 
gating is accomplished by timing delays referred to the Q restoration 
pulses. The counter system accepts a stop pulse only during a specified 
duration for each of the ten gates. These open-times must be smaller 
than the interval between Q switched pulses. If the range uncertainty is 
largejithen as well as opening the range gate s the rate at which the Q is 
restored must also be lengthened. The method of setting the range gates 
is optional, depending on the mode of operation; the minimum time delay 



can be preset for discrete times by setting the values by thumb wheel 
switches or the delay can be set in via punched paper tape at a rate of 
once or twice per second as programmed, the latter being used when 
tracking visually and firing the laser at will. For the case when only 
large uncertainty of the expected range is available, the system can. be 
operated as a single pulse system with an infinite -time range-gate; i. e. 
it is opened and remains open until a stop pulse is generated or until it 
is reset in the "no range" condition. The counter system is started and 
stopped when the signal reaches preset threshold levels - the level is set 
according to prevailing noise conditions. 

The "photographic laser", operating in the normal mode, is capable 
of about 500 joules output but can be pulsed at a faster rate with less input. 
The flash lamps are water cooled and the 330 mm by 19 mm ruby is cooled 
to liquid nitrogen temperature. The 2 ms beam is collimated to 2 mrad. 
When the output is about 200 joules the laser can be pumped at a rate of 
once every three seconds for a limited time. The PC 1000 camera used 
to photograph the returns has a 10 degree field. The approach used is to 
obtain 5 returns spaced 3 seconds apart in the same star field, swing the 
camera and obtain another set of 5 data points in the same pass. Of course^ 
range data are taken during the entire pass without interruption. Depending 
on the pass, up to three plates can be exposed. 

Several groups have succeeded in photographing Q switched laser 
pulses; it has been our experience that this present Q switched system 



is somewhat marginal. Although photographs of the return have been 
suecesful with the Q switched laser under favorable conditions, it is felt 
that with the relatively large divergence and small camera aperture 
(200 mm) that much more reliability can be achieved by using a laser 
designed specifically for the purpose. 

The two lasers, power supply and pulser for the Pockels cell, and the 
receiver are mounted in the tracking pedestal. (Modifications are also 
being made to mount a CCTV camera in the pedestal. ) The PC 1000 camera 
is located a few meters away and is operated in a non-tracking mode. The 
mount can be operated by joy stick or in a programmed mode. In the pro- 
grammed mode, five sets of azimuth and elevation angles are dialed into the 
programmer using thumb wheel switches with a resolution of 0. 001 degree. 
The pedestal automatically positions to the first set and upon receipt of a 
slew command pulse proceeds to the next position. As soon as the pedestal 
goes to position 2, the first set can be changed to the sixth set and so on. 
This allows continuous stepping of the mount for an entire pass. 

The accuracy with which the pedestal points the lasers and receiver 
is about 4 arc seconds. A print-out of the angles can also be made; this 
is especially useful when the system is operated in the joy stick mode. 

When the laser is fired a pulse is generated which is used to obtain a 
print-out of time of firing and position of the pedestal at that time. 

To operate the system it is necessary to set thumb wheel switches on 
the clock for the time a sequence is to start. (The cesium controlled clock 
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is accurate to one microsecond as compared to Loran C and several traveling 
clocks. ) Using thumb wneel switches on the "sequence controller", the 
following functions are programmed: the number of laser firings, up to 7: 

the interval between laser firings; the interval between range laser and 
photographic laser firings; the range gate (the length of time between firing 
the range laser and application of the PMT high voltage); the interval between 
Q restoration pulses; the minimum expected transit time delays; and the time 
delay to slew the pedestal to the next position-delay with respect to laser 
firing. Thenjthe width of the gate is set on the time interval measurement 
system (TIM). The width is the uncertainty of the range m microseconds 
and noise or signal pulses during this time interval will stop the counter. 

Start and stop thresholds are also set on the TIM unit. Five sets of look 
angles are programmed in the tracking pedestal control rack by means of 
thumb wheel switches. 

All data thus far collected were taken using a hand-cranked mount 
and is, therefore, limited in volume. The mount was incapable of housing 
both laser units and each has had to take its turn for testing. The most 
immediate future plans include gathering a sufficient amount of data to 
evaluate the system. Using the multipulse approach, the difference in 
range within the 300 microseconds between the first and last pulses in a 
single firing have been observed. This is another area which will be 
pursued in the near future. In order to attain more consistent accuracy 
in both range and range rate measurements, the pulse length will be 
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shortened. Some equipment modifications are also in our future plans; 
they include: (1) modification of the tracking pedestal programming 

method. It is intended to program the mount with punched paper tape 
with provisions to manually offset from the program in order to "update" 
at the last minute, and, possibly, provide for a search mode in which the 
pointing will deviate around the program in an orderly fashion. (2) 
Cooling the range laser is in the immediate future. (3) Provide 
temperature control for the narrow band filter. (4) In addition, complete 
many small modifications to improve the system. 
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A PHOTOGRAPHIC TECHNIQUE FOR 
IMPROVED LASER-RANGING ACCURACY 


C. G. Lehr, M. R. Pearlman, and J. L. Scott 
ABSTRACT 

The instrumental errors of the ruby-laser systems of Smithsonian 
Astrophysical Observatory (SAO) can be reduced from 1 to 2 m to 20 to 50 cm 
by applying corrections to the readings of the time-interval counters. These 
corrections are obtained from oscilloscope photographs of the return pulses. 
The correction is the time between the intersection of the counter's "stop" 
threshold and the centroid of the pulse. The accuracy of the correction 
depends on how well the threshold and the centroid can be located on the 
photograph. It also depends on the accuracy with which the system's electrical 
delay can be measured. The statistical procedures that were used for accurate 
determinations of the threshold, the centroid, and the delay are described. 


The returns from the ruby-laser system at SAO's Mt. Hopkins Observa- 
tory (Lehr, Pearlman, Scott, and Wohn, 1970) exhibit a variation in signal 
strength that is typical of laser systems. The variation is due in part to the 
fact that the signal varies inversely with the fourth power of the satellite 
range. It is also due to an observed "scintillation, " or random variation in 
returns from the same satellite range. The situation is illustrated in 
Figure 1. The range measurements are affected by this variation in signal 


This work was supported in part by grant NGR 09-015-002 from the National 
Aeronautics and Space Administration. 


1 


53 


006-61 



strength when they are obtained directly from a time -interval counter that 
is started by the transmitted pulse and stopped by the received pulse. The 
resolution of the counter is 1 nsec, which corresponds to 1 5 cm. But the 
duration of the pulse is 18 nsec (full width between half-power points). Con- 
sequently, counter readings corresponding to a given range can vary signifi- 
cantly if the counter stops at different points on the pulse's leading edge. 

The point at which the counter stops is the threshold level of the counter. 

Under customary operating procedures, an attempt is made to set this 
threshold at about half the amplitude of an average pulse. The laser system 
is calibrated for such a pulse and such a setting, but errors will be intro- 
duced when the return signal differs from its average strength. Figure 2 
shows how the measured range varies for a fixed theshold setting and a 100:1 
variation in signal strength. Such variations can be expected in actual 
operation. The rms error in such range measurements is probably a little 
less than half the total 5-m variation shown in the figure. In principle, this 
error can be reduced by adjusting the threshold voltage for predictions of the 
signal strength. Such predictions are based on an inverse fourth-power 
variation with range. In practice, as Figure 1 shows, the scintillation is 
so large that such a procedure is not particularly effective. 

Even if all returning signals were of the same strength, there would 
still be errors related to the intersection of the threshold and the pulse's 
leading edge. These errors come from pulse-shape irregularities that are 
due to the random emission of electrons in the photomultiplier tube. The 
fewer the number of electrons per pulse, the more irregular the pulse becomes. 
The details of this effect have been described by Lehr, Pearlman, and Scott 
( 1970 ). Figure 3 shows sketches of typical transmitted and received pulses. 
Since the transmitted pulse retains a fixed shape and size, the time T^ 
between threshold crossing and centroid remains constant. It is the corres- 
ponding time for the received pulse that changes considerably from one 
satellite return to another. However, the pulse-to-pulse variations in T , 
the counter reading, should have signs opposite to those in T and one should 
tend to compensate for the other. The compensation is not perfect because 
there are independent errors of the order of 1 nsec in the counter readings 
and in the determinations of T£. 
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The measurement of is based on a knowledge of V ? , the effective 

threshold voltage. Since depends somewhat on the slope of the pulse ! s 

leading edge, its effective value was determined experimentally. The method 
involved finding the centroid of an irregular pulse. It also yielded an accurate 
value of the electrical delay of the system. The effective threshold, the 
centroid, and the delay are treated together in the experiment described 
below. 

The experiment consisted of 48 range measurements to a ground-based 
reflector at the Mt. Hopkins Observatory. The reflector was 776. 329 m from 
the intersection of the laser system's azimuth and altitude axes. Twelve 
returns were obtained at each of four settings of the threshold voltage. 
Neutral-density filters reduced the return from the target to a strength of 
30 electrons, the level of a typical satellite return. The calculated two-way 
travel time over the 776. 329~m distance is 5180.4 nsec. This is the time 
between the peaks of the transmitted and received pulses. The system delay 
is the value that the laser system gives for K in Figure 3, less 5180.4 nsec. 
This delay is subtracted from satellite ranges that are measured to the 
centroid of the returned pulse. 

Let us consider ground- reflector returns that were obtained for a given 
setting of the threshold voltage. Since the effective threshold voltage V- 
is unknown, its value has to be determined. The determination is performed 
by as suming a number of values of threshold voltage. If we let V designate 
these assumed values, we obtain results corresponding to those sketched in 
Figure 4. For simplicity, this figure is drawn for only 5 laser returns, but 
12 returns were obtained in practice. Seven assumed values of V are shown 
in the figure. Thus, for each of the five returns we have seven values of T„, 
and one value of T , the counter reading. Figure 4 shows 35 values of 
T c + T^ plotted against the 7 assumed values of V. According to Figure 3 , 

T^ + T^ equals K when V = V^. When V = V-,, T^ and T^ should have the 
greatest correlation, so the variation in T + T_ should be smallest. This 
fact is used to find the unknown value of V-,. It is that value of V that gives 
a minimum dispersion of T + T^. In other words, it is the abscissa in 
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Figure 4 where the curve "necks down. " The corresponding ordinate is K, 
the quantity from which the system delay is determined. Figure 4 is just an 
illustration. In practice, is found by computing the correlation coefficient 
of the measured values of T c and and looking for the point where the 
correlation, which is negative, has its greatest absolute value. 

The measurement of T^ values from the photographs of the return pulses 
requires that the centroid of the irregular pulse be located. Figure 5 shows 
a triangular-shaped overlay used for this purpose. The sides of the triangle 
are separated at their midpoints by 18 nsec, the duration of the transmitted 
pulse. The absolute values of the slopes of the two sides are equal to each 
other and to the maximum slope of the pulse's leading edge. The overlay is 
moved horizontally until the sides of the triangle and the sides of the pulse 
are in best alignment, as judged by eye. Other methods of centroid deter- 
mination were tried, but this one, although simple, turned out to be quite 
accurate . 

Figure 6 shows experimental results for a threshold setting of 1 . 5 v. 

The correlation coefficient of T c and T^ and also the average of T^ + are 

plotted against the assumed voltage V. The statistical sample consisted of 

12 returns. The effective voltage came out to be 1.0 v, and the corresponding 

K value was 5268.4 nsec. The standard deviation of T +T 0 at V_ was 

c 2 2 

2. 2 nsec, which when multiplied by 15 cm/nsec gives a range error for a 

1.5-v threshold setting of 33 cm. The same procedure was repeated for 

threshold settings of 0. 5, 1.0, and 2. 0 v. In each case, V K, and the 

standard deviation of T + T-, were determined. The four values of K were 

c 2 

within an interval of 0. 3 nsec. Their average, 5268.4 nsec, was then used 
to correct the values slightly. The corrected value, rather than the one 
corresponding to the actual minimum of the correlation coefficient, is shown 

in F igure 6 . 

Figure 7 shows how the error in + T^ varies with threshold setting. 

This error lies between about 1 and 4 nsec or about 20 to 50 cm. It provides 
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an estimate of the instrumental accuracy that can be obtained from photographed 
pulses. As might be expected, it is largest at low threshold values, where the 
slope of the pulse's leading edge is smallest. 

This method was used to determine several satellite ranges from photo- 
graphed returns. In Table 1, the values are compared with the corresponding 
ones obtained from the counter readings alone. The threshold was set at 
1 . 0 v. The standard delay, which is used when photographs are not available, 
corresponds to a typical point on the leading edge of an average pulse. Its 
value was 60 nsec for the returns under consideration. The system delay 
that came from the tests described above is 88.0 nsec. Since it corresponds 
to the pulse's centroid, rather than to a point on its leading edge, it is larger 
than the standard delay. Table 1 indicates that centroid detection can 
appreciably reduce the error in the range measurement. 
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Table 1. Comparative range measurements on satellite 6503201. 
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Figure la. Laser returns. 
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Figure lb. Laser returns. 
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Figure 3. Transmitted and received pulses. 
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CORRELATION COEFFICIENT 




Figure 6. Threshold setting 1 . 5 v. 
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Figure 7. Range error vs. threshold setting. 


66 


14 




SOME COMMENTS ON THE PRESENT AND FUTURE VALUE 


OF GEOMETRIC SATELLITE GEODESY 
K. Lambeck 


T o be presented at the Geos 2 Review Meeting 
Goddard Space Flight Center, Greenbelt, Maryland 


June 1970 


Smithsonian Institution 
Astrophysical Observatory 
Cambridge, Massachusetts 02138 

67/68 


006-63 



SOME COMMENTS ON THE PRESENT AND FUTURE VALUE 
OF GEOMETRIC SATELLITE GEODESY 

K. Lambeck 

1. PRESENT STATUS OF GEOMETRIC SATELLITE GEODESY 

Geometric satellite geodesy using simultaneous optical direction obser- 
vations is a very powerful method for obtaining accurate information on 
relative station positions in the form of directions between stations. This 
method has the advantage over dynamic techniques in that it is independent 
of any orbital theory; consequently, a minimum of assumptions have to be 
made when relating physical observations to the mathematical model. On 
the other hand, simultaneous directions do not give information on the scale 
or the translation of the system, and unless the station positions are geo- 
metrically well distributed, the propagation of the uncertainties in the 
station coordinates increases rapidly if they are not constrained by some 
form of scale control. This is seen in Figure 1, which graphically represents 
the results obtained from the geometric satellite solution that is part of the 
1969 Smithsonian Standard Earth (Gaposchkin and Lambeck, 1970). The 
origin has been fixed at station 9010 (Florida), and the scale has been fixed 
by the distance 9001 (New Mexico) to 9010. Despite the large uncertainties 
in some coordinates, such as stations 9005 (Japan) and 9002 (South Africa), 
the accuracy of directions between stations is still very high. This is illus- 
trated in Table 1 for some selected lines where directions have not been 
directly observed but have been computed from the overall geometric solu- 
tion. Table 2 shows some directional accuracies for directly observed lines. 


This work was supported in part by grant NGR 09-015-002 from the National 
Aeronautics and Space Administration. 



These accuracies approach 1 prad for lines for which many data are avail- 
able. F or a 2000-km line this represents about 2 m. For numerous lines 
SAO has large numbers of unreduced data that could provide, when reduced, 
directional accuracies of this order. Comparisons with these attained 
accuracies and with the theoretically expected values have indicated that the 
influence of systematic or model errors in the data has been reduced very 
considerably. One important feature of the SAO geometric solution that 
helped reduce some of the temporal systematic errors — for example, star 
catalog errors and refraction — is that the observing campaign was carried 
out over a time interval of at least 2 years, and in many cases 5 or 6 
years. Thus, SAO has built up a data set that will provide a very valuable 
basis for future observing programs. 

The importance of the geometric solution in establishing a geocentric 
reference system is illustrated by the results of the 1969 Smithsonian 
Standard Earth. In this new solution a comparison was made between the 
station positions determined by the geometric and the dynamic methods 
before a final combination was made. These comparisons enabled us to 
verify the accuracy estimates made for the individual solutions and to adopt 
an optimum weighting procedure in the final combination. Based on the 
results of these comparisons, we feel justified in saying that the position 
accuracy of the fundamental Baker -Nunn network is better than 10 m. The 
importance of the geometric solution is further emphasized when we compare 
the results of the new Standard Earth with those of the 1966 Standard Earth 
(Lundquist and Yeis, 1966). Table 3 gives the coordinate differences. For 
the stations where it was possible to make comparisons and combinations in 
both solutions, the difference between the two Standard Earth solutions lies 
within the accuracy claims made for the respective solutions. 

For stations where the 1966 geometric solution was poor (the Z compo- 
nents of stations 9007 and 9011) or nonexistent (all coordinate components 
of stations 9002 and 9003), the difference between the two Standard Earth 
results is larger. For stations 9007, 9011, and to a lesser degree, 9002, 
considerably more data were used in the new solution, and comparisons 
between the geometric and the dynamic solutions support the accuracy claims 
for the new solution. 
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Some comparisons are given in Figure 2, which shows projections of the 
station-station vector — as determined from the geometric and the dynamic 
solutions — onto a plane normal to this vector. Both vectors have been con- 
strained to pass through the same point at one of the stations so that the dif- 
ference in the diagrams between the geometric and the dynamic solutions 
reflects the uncertainties in the coordinates of both stations. 

In addition to contributing to the coordinates of the global system, the 
geometric solution assists in separating two important and correlated factors 
in the dynamic solution: namely, the difference between the computed and 
the observed satellite positions is a function of uncertainties in the station 
coordinates and certain terms in the earth's gravity field. These differences 
are analyzed in order to compute corrections to the coordinates and to the 
gravity field. From combination of the geometric coordinate information 
with the dynamic results, an improved separation of these two factors is 
possible, leading to a gravity field that gives better comparisons against 
surface gravity than does the gravity field determined solely from the 
dynamic solution. 

2. FUTURE REQUIREMENTS FOR GEOMETRIC SATELLITE GEODESY 

The orbital theories developed at SAO include all perturbations with 
magnitudes greater than 0. 5 m. Thus, laser data of 1-m accuracy will 
enable station coordinates to be determined with an accuracy also of about 
1 m. 


At this level of accuracy and sophistication of orbital theory, the appli- 
cation of geometric satellite geodesy of comparable accuracy is essential 
for both comparison and combination purposes in order to detect and elimi- 
nate any uncertainties in the results and to separate the orbital perturbations 
caused by uncertainties in the station positions and in the earth's gravity 
field. 
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The geometric solution to date has been almost entirely based on direc- 
tion observations, but some initial work on the combination of laser range 
data with optical direction data has been carried out in France and at SAO. 

The results have already indicated the importance of this method for obtain- 
ing station positional information that is of an accuracy comparable to 
the expected results from dynamic satellite geodesy. 

The simplest approach is to observe the satellite simultaneously with a 
combination of lasers and cameras from two stations: lasers and cameras 
at both stations, cameras at both stations but a laser at only one, or lasers 
at both stations and a camera at only one. An optimization of the many 
variables that enter into solutions of this kind has already been discussed. 

The feasibility of such a program has been shown in terms of available satel- 
lites, instrumentation, and observing time span (Lambeck, 1 968). But one 
of the problems of combining laser results with optical data is their differ- 
ence in accuracies. Laser ranging now gives accuracies of about 1 to 2 m, 

or better than 1 in 10^, whereas directions can be obtained with an accuracy 

5 

of about 1" to 2", or 1 in 1 0 . Another problem is to obtain exactly simul- 
taneous observations. The partial answer to the first problem is related 
to the second - namely, the use of curve -fitting methods, interpolating for 
a fictitious instant. This reduces random errors in the data and ensures 
strict simultaneity. For the Baker-Nunn camera, however, the accuracy 
of such a fictitious observation is currently about 1'.'5 and can be improved 
to about 1!'0 with minor modifications in the reduction process (Lambeck, 

1 969 ). This still represents an accuracy of only about 5 in 1 0 , as compared 

6 

to 1 in 1 0 for laser range, data. 

But there is still another possibility for making optimum use of the laser 
data. First, we must employ directional data between the stations; these 
data, determined by simultaneous direction observations collected over the 
years, can reach accuracies approaching 1 in 1 0 . This information can 
be introduced into the adjustment of simultaneous laser optical data as 
quasi-observed quantities. Second, we must observe the satellites in 
favorable positions such that the variances in the direction observations have 
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a minimum influence on the accuracy of the desired distance of the station- 
station vector. 


For example, if the satellite P g is observed by cameras from stations 

P and P. and by a laser from station P , the distance L of the station— station 
a b ’ a 

vector is related to the observed quantities by 


sin (9 + 0 , ) 

t 2t b 

L = r : — 5 , 

a sm 0 , 
b 

where 9^ is the angle subtended at point P^ by the other two points in the 

triangle P P, P , and r is the observed station— satellite range from P . 

a d s a a 

Applying the law of combination of variances gives 




2 

sin 9 


+ cot (9 + 9, ) cr Q 

a b 0 


sin^ 9, 


sir, 2 ( e a + e b ) 



; 


and the problem becomes one of finding, for given variances of the observed 

range and the observed angles 9 , 9, , the geometrical configuration required 
^ a. d 

for (<r ) to be minimum. Several possible examples are discussed in the 
next section. 


Geometric satellite geodesy using laser range observations may also 
take the form of trilate ration — simultaneous range observations from at 
least four stations. The advantage of this method is that no directional data 
are required, so the full accuracy potential of the laser is used; with laser 
data accurate to about 2 0 cm, relative station positions can be determined 
with an accuracy of perhaps 10 cm. 
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As an example, given four stations situated on the corners of a 1000-km 
square and a satellite of height H, with n observations uniformly distributed 
through the common visibility area of the four stations, the relative accuracy 
of the station positions can be expressed by 

o- 2 ( L = 1 000 km 

— m , < 500<H<15 00 


for a maximum zenith distance of observation of 6 0°. For z = 75° 

max 


L 


260 


L 

H 


0. 55 


L \ tr -■ (L. = 1000 km 

77 ) — - m , <500<H<1500 km 

H / n iz = 75° . 

f max 

For cr = cr =20 cm, about 120 observations are required in the former 

X I* 

case, and a.bout 25 in the latter. 

For global geodesy, however, the method does require a dense distri- 
bution of stations, particularly with the satellites currently in orbit, and 
its application appears to be limited to areas where continental drift or 
large-scale crustal motions are suspected. Laser- reflector- carrying 
satellites with altitudes of about 3500 to 4000 km are very desirable for this 
purpose . 

More important, trilate ration does not give any absolute directional 
information, so it cannot contribute to kinematic studies of polar motion or 

earth rotation. 

Several combinations of laser and direction observations are possible, 
each of which has some advantage, depending on the situation in which they 
are employed. In the following section three instances are treated. 
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A. Range observations from both stations, directions from one station, 
and the direction between the two stations are known with a high degree of 
accuracy. 

B. Range and directions from both stations are known, but there is no 
a priori knowledge of the direction of the stationr-station vector. 

C. Range and directions from both stations, as well as the direction of 
the station— station vector, are known. 

3. AN OPTIMIZATION OF PARAMETERS FOR SIMULTANEOUS 
RANGE AND DIRECTION OBSERVATIONS 


Case A 


The observed quantities are the distances r^ and r^ and the direction 

P P [6 , (a - 0) 1 . The direction P P, has been determined a priori 

a s as as J a b 

from, for example, simultaneous direction observations. The angle 0 is 

St 

given by 


cos 0 

a 


i . i +m,m +n.n 
ab as ab as ab as 


with 


I st \ 

"cos 6 cos (a - 0)” 


cos 6 sin (a - 0) 

\n / 

sin 6 


{ 1 a) 


(lb) 


The distance L can be expressed as 

L = r cos 0 +r, cos 0, , (2a) 

ct sl b fo 

with 


0. = sin 



sin 0 


f2b) 
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Applying the law of propagation of variances to equations (2) yields the 

accuracy S with which L can be determined as a function of the accuracy of 

the observed quantities. If we assume that the two ranges are measured 

2 2 2 

with equal accuracy, i. e. , cr = a- = cr , then 

r r , r 

a b 


_2 

S = 


cos (0 + 0, ) + 1 

cl D 

_ 

cos 0, 


_2 . 2 si " 2 ' 6 a + V 2 

^r + r b 2 ^0 

cos 0. a 


Differentiating equation (la) yields 


sin 0 d0 = 
a a 


(lmn) as ( dm ) + {imn) ab ( dm ) 

Vdn/ ab Vdn 4s 


( 3 ) 


and if we assume that the accuracy of the direction P P can be written as 

a s 


2/10 


as \ 0 1 


and the accuracy with which the direction P P. is known as 
7 a b 


2 / 1 0 


ab 0 1 


then propagating variances through equation (3) gives 


2 1 , 2 , 2 , 

"0 “ (lr as + ^ab • 

a sin 0 


Hence , 


cos (0 + 0. ) + 1 

v a b 

2 ft 

cos 0, 


2 , r l sin2 (0 a + e b' , 2 , 2 , 

T r + TZ 7T— "as + "ab* • (4a) 

cos 0, sm 0 
b a 
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or, alternatively. 


cos (0 + 0, ) + 1 

a b 

2 ft 

cos 0 , 


9 L 2 Sin 2 (0 + e, )(o - 2 + ( 7 2 , ) 

2 , 'a b as ab 

£ r + r 

r 2 


(4b) 


cos 


0 , 


If n sets of observations have been made, the accuracy with which the dis- 
tance L can be determined is 



S will be a minimum when 0^ = 0° or 180° and will be infinity when 

0, = 90° - 

b 


The curvature of the earth, the satellite height, and the maximum 
zenith distance at which observations can be made impose limits on the pos- 
sible values for 0 and 0, . To express 0 and 0 as functions of zenith 

cl 0 3. D 

distance z, satellite height H, subsatellite distance r), and azimuth (3 of the 

3 . 

subsatellite point relative to P^P^, we use the following relations (see 
Figure 3), where R is the radius of the earth: 

r a r b 

sin p a = sin z a • ; sin ^ = sin z b • ^ , 


ib _ L 

sm 2 = 2R * 


cos T), - cos 4* cos r| 

0 b a 

COS p = : i : 

a sm t sin i] 


; cos (3 


COS T) - COS Ip cos 9, 

a b 

sin Lp sin p. 


ib dj 

cos 0 = cos p sin z cos — - cos z sin — , 

a a a 2 a 2 


ib .ip 

cos 0, = cos p, sm z, cos v- - cos z. sm -~ 
b b b 2 b 2 
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2 

The function S is evaluated as a function of r| and (3 in Figure 4 for 

2 2 2 1 /2 

some typical values of L and H for t r =2 m, cr = (or + cr ) = l'.'O, and 

r s 3 , s 3 . b 


max 


75° and z 


max 


6 0 ‘ 


For a given maximum zenith distance at which observations can be made 
and for = 0, cos 0^ is largely independent of the distance L (or angular dis- 
tance 4*)- This is verified when the minimum value of is computed for 

satellites with different heights and for various values of L. Figure 5 gives 

Z 

these minimum values of S as well as the number of observations required 

to determine the distance L with an accuracy equal to cr^ if the observations 

are made when the satellite is near the optimum positions. The results are 

given for 1) z = 75° cr =2.0 m, cr^ = 1!'0; 2) z = 60° , cr = 2. 0 m, 
max ’ r S max r 

o-“ = 1 ! 1 0 ; 3) z = 75°, o- = 2. 0 m, cr = 2!'Q; and 4) z = 75°, <r = 1. 0 m, 
s max r s max r 

0 -2 = i"o. Clearly, it becomes desirable to observe as close to the horizon 

S 

as possible at all times and to measure the angle 9 with as great an accuracy 
as possible. An improvement in the accuracy of the range observations does 
not significantly reduce the number of observations required. 


With Baker-Nunn cameras the accuracy of an interpolated direction 

observation is currently about 1 1 ' 5 but can be improved to about l'.'O. Then, 

2 2 

if the direction is known a priori with high accuracy, & s ~ 1.0 arcsec . 

If the direction has not been determined from simultaneous observations 

but can be computed from the coordinates derived from a dynamic determina- 

2 2 

tion, then cr is of the order of about 1" and cr = 2. 0 arcsec . The dynamic 
J as s ’ 

and geometric solutions for the length L would, however, no longer be 

independent. 

If there are cameras at both stations, the satellite can be observed 
simultaneously with both lasers and cameras; but the optimum configuration 
will be different, since 0^ is now an observed quantity. This is discussed 
further below. Alternatively, we can observe the satellite with one camera 
and the two lasers when the optimum configuration occurs for this combina- 
tion, or we can observe the satellite with the two lasers and the other camera 
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when the appropriate configuration occurs. The advantage of this alternative 
approach is that the satellite need be visible only from the camera site if the 
laser at the other site is powerful enough to obtain daylight returns. 

The requirement that the satellite be observed in a very limited part of 
the sky is rather stringent and could result in a long observing program. 

But this should not be a real problem, since only a small number of obser- 
vations is required, the laser at can be operated in daylight if it is power- 
ful enough, and several retroreflector-carrying satellites are available. 

This point is illustrated further by examples in Section 4. 

Case B 


Let us consider now the case where the satellite is observed simultane- 
ously with both cameras and lasers from the two sites and where no precise 
information is available on the relative stations positions. Then, 


where 


r , 2 , 2 _ „ .1 / 2 

L = (r + r, - 2r r, cos 0 ) 
a b a b s 


cosS - 1 i. + m m. + n n, 

s as bs as bs as bs 


With 


2 1 


• 2 ft 

sin 0 


, 2 , 2 . 
(a as + °bs> 


2cr 


• 2 a 
sm 9 


if 


2 2 , 2 , 

%s = ^bs = • 


(5) 


1 1 
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propagating variances through equation (5) gives 



For given L and H, S is a minimum for a satellite position midway 

between and P, (see Figure 6), and the larger the distance L, the smaller 

2 a ° 2 
will be S . The relation between L, H, and S . is given in Figure 7 for 

or = 2. 0 m and cr^ = l'.'O. We see that for this case, . is very dependent 

on both L and H. On the other hand, the results in Figure 7 are independent 

of z , at least in the range of 6 0° to 75°, since when the satellite is 

max ° 

midway between the two stations, it is generally at a zenith distance of less 

than 6 0° . 


For the Case B combination of observations, the satellite must always 
be visible from both stations and daylight observations are not possible. 


Case C 


If, in addition to both directions and ranges, the direction of the station- 
station vector is known, we obtain two estimates for the length L: 

r sin (0 + 0, ) 

T _ __a a b 

J 1 sin 0^ ' 


L 2 = 


r, sin (0 + 0, ) 

b a b 

sin 0 


and an adjustment must be made to obtain a unique and most probable value 
for L. The accuracy estimate for this adjusted length is given by 
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s 2 = 



4 4 

_ , sin 0 + sin 0, 

cr + [l + cos (0 + 0, )] 

r L ' a b ,J 


. 4 „ .4 

sin 0 sm 
a 


ST 


+ 


T 2 

4 } 

L 

(T 1 


s l 

• 2 o 
sm 0 

sin 2 0,\ 

a 

b ) 



[cos (0 + 0, ) sin 0 + sin 0,1 

1 a b a cr_ 

■ 4 a 
sm o 


+ 


[cos (0 + 0, ) sin 0, + sin 0 ] 

L v a W b £r_ 

• 4 ft 

sm 0, 


ST 


sin (0 + 0, ) 

a b ' 


with 


2 2 2 
cr r = <r r = cr 

a b 


2 2 2 2 2 

cr = cr + cr , = cr. + <r, 
s a ab b ab 


,2 . 


In Figure 8, S is evaluated for L = 1500 km, H = 1000 km, by use of 


or = 2. 0 m and cr = 1 ! 1 0. 
r s 


.2 . 


The function S is a minimum for a given L and H when either 0 or 

2L 

is as large as possible, which occurs when (0^ + 0^) is a maximum — i. e. , 
when the angle subtended at the satellite by the two stations is a minimum. 


The minimum values of S as a function of L and H are given in Figure 9 


for z =75 
max 


For z = 60°, the relationship becomes less dependent 

ii 1 cL X 2 

on H. The average values of as a function of L are indicated by the 
broken line in Figure 9- 
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4. SOME EXAMPLES 


As part of the ISAGEX observing program scheduled for 1970-71, 
it is intended to obtain simultaneous range and direction observations between 

several stations: 

8015 Haute Provence (laser) 

9004 San Fernando (laser and Baker-Nunn) 

9029 Natal (laser and Baker-Nunn — daylight laser observations 

will be possible) 

9091 Athens (laser and Baker-Nunn) 

9041 Dakar (laser and Baker-Nunn) 


Determination of the following distances will be attempted: 


8015-9004 1306 km 
8015— 9091 1656 km 
9004-9091 2648 km 
9004-9041 2600 km 
9029—9041 3000 km 


From the 196 9 Smithsonian Standard Earth results and from the Centre 

t 

N ational d 1 Etudes Spatiale s RCP- 1 33 results, the directions of the triangle 
801 5 — 9004 — 9091 have been accurately determined by the geometric method. 
For lines 8015—9004 and 8015—9091, Case A is therefore appropriate, while 
for the line 9004—9091, Case C is the most suitable. For the lines 9004— 
Dakar and 9029—9041, Case B has to be used, since no directional informa- 
tion on these vectors is currently available. Figure 10 gives the number of 
observations necessary to determine the statiormstation vector with an 
accuracy of 2 m, assuming the optimum conditions are satisfied. Curves 
for three satellites are given. These and others can be constructed from 
the general cases given in Figures 5, 7, and 9- If o' = 1 1 1 0, cr =2.0 m, 

S I* 

and z = 75°, only 14 observations are required to determine the dis- 

max 

tance 8015—9004 with an accuracy of about 2. 0 m, and 18 observations are 
required for the distance 8015—9091. For these lines any satellite at any 
height can be observed with equal effectiveness. For the line 9004—9091 the 
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number of observations required varies between about 10 and 16, depending 
on the satellite height. For the remaining two lines the situation is less 
favorable. In both instances the satellites must be observed when they are 
near perigee, as height now plays a very critical role. For example, using 
Geos 1 for the line 9004—9041, about 94 observations are required if the 
satellite is observed near apogee, whereas only 16 observations are required 
when the satellite is near perigee. A similar number is required for 
Geos 2 observations at perigee. 

The long line Natal— Dakar will present some problems, since observa- 
tions near perigee will not always be possible because of the earth's curva- 
ture. F or a height of about 1200 km, at least 20 observations are required. 

For both these lines the situation will be considerably improved if the 
station— station direction is also known (Case C). 

5. CONCLUSIONS 

For a successful program of simultaneous direction and range obser- 
vations, the following points must be considered: 

1. The satellite— station configuration is extremely critical in Cases A 
and B if the number of observations are to be kept to a minimum. The 
conditions are such that a particular satellite satisfies them only very 
occasionally, but with the six retroreflector-carrying satellites now avail- 
able, sufficient data can be collected in an observational period of about 

6 months. 

2. It is important that the direction of the station— station vector be 
independently observed by simultaneous direction observations with an 
accuracy of about 0!'5 or better. This can be done by use of more favorable 
satellites and the optimizations discussed by Lambeck (1969). 
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3. The direction observations should be accurate to about 1!'0 or better. 
With care, this can be achieved for synthetic Baker-Nunn observations, 
particularly as the lines discussed here are mostly in the Northern Hemis- 
phere and stars of far south declinations are not required in the film- reduc- 
tion process. 

4. .Laser range accuracies of about 1 or 2 m are adequate for this 
observational program. Improved range accuracies do not significantly 
reduce the number of observations required. 

5. If the interstation distances are to be determined with accuracies 
better than about a meter, range observations only in the form of spatial 

tr ilateration should be used. High-altitude retr or eflector- carrying satellites 

will be required for this. 

6. For the 1 -m range accuracy and for low laser-firing repetition rates, 
interpolation methods for obtaining synthetic simultaneous observations are 
adequate if the direction observations are taken on either side of the range 
measurement and if the interpolation is carried out at the instant at which 
the latter observation is made. 

7. For spatial trilateration either very rapid repetition rates are 
required or the lasers must be synchronized so that exactly simultaneous 
observations are possible. For 20-cm accuracies, this means that the 
absolute time at each station must be known to within about 2 0 psec and that 
the laser can be fired on command. Lasers employing Pockel cells for 
triggering the firing are appropriate for this. 
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Table 2. Direction cosines of station— station vectors. 


Station 

L 

M 

N 

2 

*6 

2 

a 

a 

2 

900|9007 

*5533024 

-.1013374 

-*8267933 

.08 

.04 

e * © 02 

90019009 

,8673538 

-.1488313 

-.4749175 

.04 

.05 

-.0? 

90019010 

,9654354 

-.1669486 

-.2001564 

.05 

.06 

-.03 

90019012 

-,7952943 

.5590340 

-.2344950 

.22 

.26 

.16 

90019113 

-.8398600 

.4984305 

.2149468 

.55 

.44 

.18 

90019114 

.1092630 

.6856665 

.7196687 

.48 

.16 

-.02 

90019117 

-.7167568 

.6500002 

-.2525064 

.22 

.29 

.15 

90029008 

-.2634801 

.2647682 

.9276185 

.26 

.86 

— .41 

90029028 

-.0386285 

.3164744 

.9478142 

.65 

.47 

.17 

90049006 

-.5590305 

,8242178 

-.0902773 

.14 

.17 

-.10 

9004900 S 

-.3267921 

.9374792 

-.1197483 

.20 

.13 

-.11 

90049009 

-,4414247 

-.8138821 

© 3778096 

.31 

.37 

,25 

90049010 

-.6274848 

-.7660084 

-.1351580 

.29 

.40 

© 2 2 

90049028 

-.0379168 

.8490277 

-. 526-9860 

.38 

.38 

-. 06 - 

90049029 

.0149787 

-.5736286 

-.8189785 

1.09 

.38 

.30 

90049091 

-.1927372 

.9797634 

.0540006 

.31 

.31 

-.15 

90049115 

-.6890483 

.3985934 

.6052568 

.72 

.51 

© 2 1 

90059006 

.9152361 

.3880019 

-.1086161 

.52 

.42 

.38 

90059012 

-.2473509 

-.9394576 

-.2371432 

1.06 

.27 

1 © i 9 

90059117 

-.3907736 

-.8491932 

-.3551997 

2.04 

.42 

1.62 

90069008 

.9110432 

-.4121827 

.0102782 

© 6 2 

.28 

.23 

90069028 

,8289764 

-.3212882 

-.4577905 

.68 

.44 

.16 

90069115 

.3606918 

-.8366851 

.4121402 

.23 

.37 

,15 

90079009 

.0984437 

-.0040855 

.9951342 

.05 

.10 

.03 

90079010 

-.2021848 

.0424032 

.9784290 

.05 

.06 

® 02 

90079011 

.1850049 

.4871412 

-.8535026 

,22 

.15 

©08 

90079029 

,7997417 

.5301435 

.2817112 

.34 

.74 

-.03 

90079091 

-.0766857 

.5210859 

-.8500522 

.65 

.58 

.17 

90069028 

.5673330 

-.1630395 

-.8071873 

1.47 

.09 

.28 

90089115 

-.0568149 

-.8471896 

.5282441 

.46 

.20 

.14 

90099010 

-.6310609 

,1066281 

.7683701 

.10 

.17 

.04 

9OO99OI 1 

,0060328 

.1892182 

-.9819165 

.08 

.03 

.02 

90099029 

.7072600 

.5213013 

-.4775230 

.75 

.75 

© 0 1 

90109114 

-.5807344 

.5531073 

.5973440 

.16 

.14 

.06 

901 19029 

.6980520 

.3028520 

.6488482 

.50 

.81 

-.10 

90129114 

.8019840 

-.2028475 

.5618492 

.46 

.40 

« © 00 

90129117 

-.3703243 

.8841374 

-.2848877 

1.70 

.43 

© 9 ^ 

90289091 

-.0872703 

-.5447064 

.8340736 

1.22 

.90 

-.03 

90299031 

-.6643695 

-.0872138 

-.7422984 

.43 

.73 

-.04 

91139114 

.5223308 

.5101516 

.6833124 

.55 

,45 

-.12 

91149115 

.7328034 

.6782404 

.0546731 

1.22 

.74 

.20 

80159004 

. 4 Q 36878 

-.7757363 

-.4850458 

.51 

.32 

.18 

80159066 

-.6962356 

.3087623 

.6480138 

8.83 

.68 

1.29 

80159074 

-.7231334 

.4996525 

.4768914 

1.31 

.42 

-.53 

80159080 

-.6121440 

-.5512713 

.5669035 

1.42 

.65 

.35 

90049066 

-.4796719 

.6955520 

.5349041 

.34 

.18 

-.08 

90049074 

-.6073169 

.6246964 

.4908366 

.44 

.14 

© 1 7 

90049080 

-.6703400 

.2377845 

.7029245 

.45 

.19 

-.01 

90669074 

-.7225951 

.5377779 

.4343400 

2.01 

.66 

-.75 

90669080 

-.4578721 

-.7820683 

.4227555 

1.66 

.48 

.49 

90749091 

.6757205 

.2958945 

-.6751653 

1.91 

.55 

@21 

90021043 

-.4604214 

.8202075 

.3395169 

1.16 

0.89 

1.09 
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Table 3. Differences in the coordinates as determined in the 1966 
and 1969 Standard Earth solutions. The large discrep- 
ancies occurred at stations where no comparisons or 
combinations of independent solutions were possible in 
the 1966 solution. 


Station 

X 

(m) 

Y 

(m) 

Z 

(m) 

9001 

+ 1.7 

-0.6 

+ 1.0 

9002 

-0.7 

+26. 0 

+32. 1 

9003 

-26. 0 

-13. 8 

+26. 5 

9004 

- 4. 7 

+ 3.7 

- 7. 1 

9005 

+ 4. 3 

+ 13. 4 

-11.3 

9006 

- 2. 2 

+ 3.3 

+ 8.6 

9007 

+ 5.6 

- 3.4 

+ 27. 8 

9008 

+ 11.2 

- 9.0 

- 4. 2 

9009 

+ 8.6 

- 3. 6 

- 4. 1 

9010 

+ 9.2 

- 9.2 

- 2. 5 

9011 

+ 14. 2 

- 4.0 

+30. 9 

9012 

+ 1.8 

- 6.8 

+ 0.7 

9114 

+ 1.4 

+ 8.4 

- 8.0 

9115 

+ 4.9 

+ 14. 6 

+ 5.6 

9117 

+ 3. 1 

+ 13. 9 

+ 4. 5 
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igure 1, Error ellipses for the Baker-Nunn camera coordinates derived from the geometric 
solution. The origin is station 90 1 0, and the scale is the distance 9001—901 0. The 
number in parentheses at each station refers to the standard deviation of the height 
dete rmination, 
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Figure 2. Representation of some comparisons for directions between sta- 
tions derived from the geometric solution A , the dynamic 
solution □ , and the combination solution □ . 


90 


22 




23 





1000 2000 3000 


2 

i? igure 5. Evaluation of S m ^ n for Case A and the minimum number of obser- 
vations required (N) to determine the distance between the stations 
with an accuracy of 2 m. 
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Evaluation of S for L = 1500 km, H = 1000 km for the case of 
range and direction observations from P a and Pj-, and with pro 
cise direction information for the vector P & Pb (Case C). The 
limits of the common coverage area are shown for z = 60° 
and z = 75 
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igure 10c. Number of observations N required for satellite 6800801. 
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GODDARD MOBILE LASER (MOBLAS} SYSTEM DESCRIPTION 

by Don A. Pretno 

Introduction 

The Goddard Mobile Laser Tracking System (KCBXAS) is a pulsed ruby 
laser radar designed primarily to track satellites equipped with optical 
corner reflectors. The system began operations in late 1968 at Goddard 
along side the Goddard Experimental Laser Tracker with which it was 
int er compared . It was during this initial series of tests that the 
MOBLAS first demonstrated the ability to track satellites in daylight. 
Since then, the tracker has been moved to Carnarvon, Australia; Mt . 
Hopkins, Arizona; and back to Goddard for collocation tests and participa- 
tion in various world-wide concentrated tracking efforts. Presently, 
the MOBLAS is in New York state for three months in support of the 
Preliminary Polar Motion Experiment (PME) . 

Figure 1 shows the MOBLAS installation at Carnarvon. The Tracking 
Pedestal is visible at the far right and the Control Van can be seen 
just to the left of a small building near the Tracking Pedestal. Not 
shown is a Power Van containing diesel generators for use when adequate 
power is not available. 

The Tracking Pedestal can be seen in greater detail in Figure 2. 

The laser is located within the air conditioned enclosure just visible 
in the lower left corner of the picture. The laser is stationary and 
the light is directed to the telescope on the left through a coelostat, 
a portion of which is visible on the azimuth axis just above the platform 
level. The central telescope collects the reflected pulse and the 
telescope on the right is for operator use. 
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Figure 



* * 


S v stem Ci iaract eristics 

The characteristics to be described do not include improvements 
soon to be incorporated for the PPME , but are valid for all prior 
measurements. The characteristics are shown in Figure 3. 

The 0.5 joule output is obtained from a Goddard-built q-switched 
laser oscillator whose active element is a ruby rod 3/8 inch in diameter 
and nominally 6 inches long. The q- switch, a rotating prism with 
b teachable dye cell, is adjusted so that a single 18 nanosecond, 

0„5 joule pulse is emitted each time the laser is triggered. These 
parameters are measured at least once each shift by use of a suitable 
thermopile and oscilloscope and monitored at all other times by an 
oscilloscope. The laser is water cooled and has been operated up to 
10 minutes at a time at a one pulse per second rate with output characteristics 
remaining constant within 307. . The beam divergence figure is based on 
a nominal 5 milliradian total angle at the laser which is reduced to 
the 1/3 milliradian value by the 15 power transmitting optics. These 
optics are of sufficient diameter to prevent vignetting. 

The receiver aperture shown is the nominal value. The effective 
value is more nearly 15 inches. The receiving telescope is a cassegrain 
type and has a focal length of 227 inches and a field of view of 
5 milliradians . However, to reduce background noise, especially during 
daylight operations, the field of view is reduced by a field stop to 
match the transmitted beam width. Additional noise reduction is achieved 
by placing a loX filter centered at 6943$ in the optical path. This configura- 
tion has been found satisfactory for ranges out to at least 2000 km when 


operating in daylight with GEOS-B. 
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MOB LAS SYSTEM CHARACTERISTICS 


LASER ENERGY/ PULSE 
PULSE WIDTH 
REPETITION RATE 
BEAM DIVERGENCE 
RECEIVER APERTURE 
RECEIVER FIELD OF VIEW 
SPECTRAL FILTER 
RANGE RESOLUTION 
TIMING ACCURACY 


0.5 JOULE 
15-25 NANOSECONDS 
1 PPS 

1/3 MILLIRADIAN 
16 INCHES 
1/3 MILLIRADIAN 
10 % AT 6943 % 

10 NANOSECONDS (] .5 METERS) 
1 50 MICROSECONDS TO USNO 


FIGURE 3 



- 3 - 

The range resolution is limited to 10 nanoseconds by the Range 
Time Interval Counter presently being used. The accuracy of this 
counter as quoted is ± 10 nanoseconds ± time base accuracy. To date, 
the time base accuracy has been held to within 2 parts in 10 referred 
to the standard offset frequency broadcast by WWVL. 

The timing accuracy stated is conservative, insofar as the MOBLAS 
time standard is concerned. By the use of VLF measurements and flying 
clocks, it is estimated that time was known to within ± 15 microseconds 
at Carnarvon and to ± 5 microseconds elsewhere. The reason for the ± 

50 microseconds is that the laser firing time is not measured to a 
greater accuracy. 

B lock Diagram 

To understand the MOBLAS from a functional viewpoint, one might 
start with Figure 4, the Block Diagram. Most of the Figure is self- 
explanatory, the system being relatively a simple one. However, one 
should note that the Range Time Counter is activated by the leading 
edge of the start and stop pulses. The effect is that large amplitude 
variations in the received energy can and do cause the range measurements 
to appear biased. However, the magnitude of the bias is limited roughly 
to the width of the laser pulse; i„e., to less than 3 meters. In future 
versions of the system, the detection circuitry will be modified to 
overcome this source of error. 

Pointing and data handling are performed by the computer. Experience 
has shown that tracking operations, especially in daylight, are much 
improved if real time adjustments can be made to satellite predictions. 
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MOBILE LASER RANGING SYSTEM 










- 4 - 

The major source of prediction error encountered so far has been the 
uncertainty in the instantaneous value of mean anomaly. The correction 
to this problem in an operational sense has been to evaluate the predic- 
tions immediately after each pass to obtain eventually a history of 
adjustment values that can be extrapolated. 

To improve pointing accuracy , a periodic astroSomical calibration is 
performed on the instrument. The procedure is to track 40 or 50 stars 
spread throughout the visible hemisphere and to use the measurements to 
adjust the coefficients of a pointing error model to reduce in a least 
square sense the difference between observed and calculated positions. 

The rms of the residuals after the adjustment is typically from 6 to 8 arc 
seconds. The coefficients are used during a pass to adjust pointing 
angles to the correct values and afterwards to correct recorded angular 
data . 

The computer has turned out to be a powerful tool, both for the 
determination of the adjustment values and for their use. 

Data Handling 

The MOBLAS is geared to the handling of a fair amount of tracking 
data per pass. As an example, Figure 5 shows a plot of the residuals 
for two passes of GEOS-B taken at Carnarvon. The upper one contains 
457 observations, the lower, 402. One advantage is that it is easy to 
decide which points are acceptable. A second is that operational func- 
tions such as satellite acquisition and continuity of track are immeasurably 
easier to implement, thus, assuring an adequate accumulation of data. 

A third is that usually a good range of measurement geometry is obtained 


for collocation type experiments. To illustrate the accumulation of data 

110 



CARNARVON LASER RANGE RESIDUALS 






SUMMARY OF CARNARVON LASER RANGE DATA 


FIGURE 6 


j Date 

T i me 


No . o f 
Obs . 

Rms 

(Meters ) 

Date 

T i me 

: . -'V- 1 -• . ■■ 

IgjEHCBB 

No . o f 
Obs . 

Rms 

(Meters ) 

1 Feb . 6 

2320 

D 

55 

1 .06 

Apr . 5 

1 140 

N 

243 

1 .55** 

S 

956 

D 

349 

1 .29** 

6 

01 22 

D 

1 98 

1 . 24** 

7 

2340 

D 

292 

1 .20** 

6 

1 1 57 

N 

397 

1 .49 

8 

1019 

D 

122 

1 . 34** 

8 

1236 

N 

76 

1 .51 

1 5 

1040 

D 

115 

1 .76** 

9 

1106 

N 

217 

1 .31 

1 6 

1100' 

D 

146 

1 .22 

10 

0049 

D 

470 

1 .07** 

16 

2252 

D 

1 64 

1.31* 

10 

1 1 24 

N 

468 

1 . 34 

1 7 

1121 

D 

25 

0.87 

1 1 

01 09 

D 

540 

1 .07* 

17 

231 1 

D 

239 

1 .58 

1 1 

1 1 52 

N 

449 

1 .45** 

18 

2329 

D 

335 

1 .44** 

1 2 

1 38 

D 

501 

1 .00 

1 9 

1 007 

D 

252 

1 .24** 

1 4 

» 17 

D 

1 89 

1 .09 

1 9 

2350 

D 

371 

1.41** 

14 

206 

D 

1 4 

2 29** 

20 

1 026 

D 

294 

1 .52** 

1 4 

1 240 

N 

309 

1 .80 

21 

0009 

D 

224 

1 . 05** 

1 5 

35 

D 

302 

1 .22 

21 

1 044 

D 

1 59 

1.61** 

1 6 

1 1 30 

N 

268 

1 .40 

23 

2316 

D 

310 

1 .74** 

17 

1 1 32 

N 

88 

1 . 80 

25 

101 1 

D 

245 

1 .08 

21 

’0040 

D 

336 

1.17** 

25 

2354 

D 

358 

1 . 38** 

21 

1116 

N 

223 

1 . 08 

26 

1030 

D 

386 

1 . 35** 

21 

1 305 

N 

1 44 

0.97 

27 

0014 

D 

409 

1 . 09* 

22 

1 1 34 

N 

366 

1.41 

27 

1050 

D 

398 

1.41* 

23 

1 1 54 
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* - Indicates Residuals are Significantly Non-Random - 5 Percent Level. 
** - Indicates Residuals are Significantly Non-Random - 1 Percent Level. 
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Figure 6 






















- 5 - 

Figure 6 shows a summary of the data taken at Carnarvon. 

Availability of the data is certainly another factor that must be 
considered. The MOBLAS can process its own data and transmit all or 
part by teletype to the experimenter within hours after it is taken. 

This was started at Carnarvon and has been continued during ever}?- 
experiment since then. Although in earlier tests not fully corrected 
for refraction effect s, the data now is fully corrected for all known 
parameters and is as accurate as that submitted to the Data Center; i„e., 
average rms of fit to short arc of 1.2 meters. 

Future Plans 

In the near future, modifications will be incorporated to increase 
range time measurement resolution by an order of magnitude and to eliminate 
sources of bias and error such as that inherent in the present range 
measurement technique and in certain time measurements. The methods 
are known, only the implementation remains. 
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INTRODUCTION 


This study was directed towards the problem of obtaining a geodetic 
position of a transponder array in the broad ocean. The transponder 
array is used for scoring the accuracy of impact of missiles,, hence 
the location of the array must be known quite accurately. The approach 
used in the investigation was to have a ship locate itself by tracking 
a satellite whose orbit is accurately determined by land based sensors. 
Then the ship may interrogate the transponders in the target array 
and thus locate the geodetic position of the array. The problem as 
a whole becomes somewhat complex, involving the determination of the 
relative positions of the transponders and the orientation as well as 
the latitude and longitude of the array. This presentation will be 
confined to a discussion of the methods and accuracy achieved in ship 
location only. 

All data were obtained from beacon track of the GEOS-B satellite. 
The orbit of GEOS-B was determined in every case from tracking data 
obtained by AFETR pulse radars. Clearly a single observation in range, 
azimuth and elevation from the ship would serve to locate the ship. 

Such a simple non-redundant piece of data would, however, give very 
poor accuracy. In practice it is desirable to have ship radar track 
from satellite rise to satellite set. During this interval, of course, 
the ship may travel two miles in not necessarily uniform motion, and 
this motion is a very important consideration in the solution. Three 
methods of determining ship motion were investigated: (1) Direct 
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solution for average ship velocity vector in least squares solution, 
( 2 ) Use of ship navigation data (SINS) for relative ship position, 
and (3) Use of ocean bottom transponder interrogation data for 
relative ship position. 

The results of the ship location tests conducted at the AFETR 
were evaluated to an accuracy of about five feet by a network of land 
based theodolites which tracked the moving ship as the ship tracked 
the satellite. The theodolite measurements were also used to 
simulate an "ideal" transponder array to provide relative ship 
positions. The complete test design included simulation tests as 
well as tests with real data. Simulations were used to study the 
following variables: 

a. Recovery of radar range bias 

b. Maximum elevation of satellite 

c. Use of radar range-only versus azimuth and elevation with 
range data, 

d. Various portions of satellite pass. 

This abbreviated discussion will be limited to a review of some of the 
representati ve results. 



DATA COLLECTION 


The data for these tests were collected at separate periods. 

The first test was conducted July 22-23, 1969. Due to priority 
problems, no additional passes were scheduled until September 1969 . 
The second test was scheduled and seven passes were collected in 
December 1969. The latter test has no absolute reference. Table 1 
shows the revolution number and the dates of radar track. It was 
intended to collect more land based radar data for the reference 
orbit but priority prevented using more than two land based radars 
in some cases. 
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TABLE 1 
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DATA REDUCTION 


Two separate data reduction routines were evaluated for this 
study* The first method treats the data exactly like a fixed land 
based radar site. Ship's radar data are transformed to a fixed site. 
Each data point requires a relative position for the ship. The ship's 
radar data were processed for input to the NITE computer program as 
follows: 


P Tfrp 


ARED 


GEST 


Convert. Scale 


Edit. Correct 


Stabilize , correct, 
transform to fixed 
geodetic position 
(X ,Y ,Z cartesian) 


| DPGE 


Transform to 
Radar A,E,R 

The second data reduction method utilized the radar range data 
after processing through GEST, The range only ship data are input 
along with earth fixed E, F, G satellite data into the HESP computer 
program. Input information required consists of (1) satellite trajectory 
points in an earth centered, earth fixed coordinate system (2) ship radar 
range data at times coincident with the satellite trajectory coordinates, 
and (3) either relative ship position points or a priori estimates of 
ship position and velocity at epoch. 

As indicated above, if the ship does not have uniform position 
changes or velocity, relative ship position coordinates may be input 
directly and the equations of motion for the ship are by-passed and 
the adjustment process is entirely unaffected by erratic ship motion. 
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TEST RESULTS 


a, NITE Computer Progra m. The land based radar range data 
and the ship's radar range data are used to simultaneously estimate 
the orbit and to locate the position of the ship. The NITE computer 
program is capable of adjusting for a range bias. Results for this 
limited series of tests indicated no significant improvement when 
radar range bias is also adjusted. The results also indicate that a 
one revolution or two revolution fit are about equal. Since the ship 
moves during the test 9 it was not considered feasible to constrain 
the survey adjustment in the same manner as a fixed radar site. 

Results of the control tests conducted near the Florida 
mainland appear in Table 2. The second test was far out in the 
ocean and no absolute accurate reference was available. The dis- 
persion was computed about the arithmetic mean. All data were 
referenced to a fixed transponder. The results of the second test 
are summarized in Table 3. The forward motion of the ship was 
determined by ASPS interrogation for the second test (Dec 69). 
Previous survey results had determined the relative position or base 
line length between three transponders. The ASPS ranging data were 
recorded at one minute intervals and relative ship position data were 
smoothed for this computer run. 
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TABLE 2 

TWIN FALLS GEODETIC SURVEY ERROR FROM NITE COMPUTER 


Rev No, 

Single Rev 

Orbit Fit 

Two Rev 

Orbit Fit 


Delta Lat 
(feet) 

Delta Long 
(feet) 

Delta Lat 
(fe'et) 

Delta Long 
(feet) 

7156 

-28.0 

101.0 

-39.7 

53,8 

7157 

-40,0 

109.0 

-63.0 

-162.0 

7169 

-35.0 

16.5 

-2.4 

69.0 

7170 

* 

* 

156.3 

-125.4 

7708 

1.2 

-19.6 

-26.8 

-64.2 

7709 

** 566 

**1322.6 

**533.1 

**1302.5 

7715 

3.1 

-107.0 ' 

NA 

NA 

7721 

-1.7 

-32.4 

64.0 

■54.9 

7722 

-2.4 

100.0 

9,6 

128.0 

7728 

4.3 

-150.9 

NA 

NA 

7734 

96.8 

125.3 

3.4 

56.8 

7735 

-79.0 

-190.1 

-5.9 

32,2 

Mean** 

-13.5 

-4.8 

-24.1 

40.8 


51 

115 

61 

88 


* Rev 7170 was not processed as a single pass since only one 
land-based radar tracked during this pass. 

** Rev 7709 results were excluded - high PCA elevation may have 
caused a poor reference orbit. 


NITE COMPUTER SOLUTION FOR Cl 2 ARRAY 

RADAR RANGE ONLY DATA FOR GEOS-B 
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*N0TE: Mean off set from a reference transponder with the 
geodetic location established by a different system. 


ESTIMATE OF SHIP POSITION AMD VELOCITY 
WITH RADAR RANGE ONLY DATA 
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*GE0S-B orbit was determined from land-based MIPIR radar data using TRACE-66 computer program* 
**Rev 7734 is not included in evaluation due to poor ship velocity estimate. 



TABLE 5 

ESTIMATE OF SHIP POSITION 
WITH RADAR RANGE ONLY DATA 
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*GEOS«B orbit was determined from land-based MIPIRradar data using TRACE-66 computer program. 
Relative position of ship determined from theodolite. 


be HESP Computer Program , The HESP computer program requires a 
reference orbit. Input data are ships radar range and the adjusted 
parameters may be varied. We adjusted position and velocity in one 
study (see Table 4). It may be noted that the estimated velocity of 
the ship was in error but the ship's position errors are small. The 
ship location accuracy was much improved when only position was 
estimated (see Table 5). The ship's relative position was input as 
a known parameter. Relative position was based on theodolite data. 

If the ship has uniform velocity during the tracking period s it is 
not necessary to know the point-by-point position during the tracking 
period. Only a good estimate of the velocity is required. Slight 
deviations from uniform velocity will have only a minor effect. 

Data were to be processed in the HESP program for the C-12 
array using SINS velocity. A delay in data reduction has caused the 
results to be delayed. Preliminary data review indicated that ship's 
radar acquired track of GEOS-B satellite late approximately fifty 
per cent of the time. A further review of data acquisition is needed 
to resolve the problem. Previous studies have indicated that 
horizon-to-horizon track by the radar is highly desirable. The real 
ship's radar data for Air Force Eastern Test Range were less than the 
anticipated track time. This possibly could degrade the results 
below the optimum solution. 



SIMULATED DATA 

Figure 1-A depicts the location of three AFETR land based radar 
trackers and four ships. The ship location was generally chosen to 
study the effect of geometry on geodetic survey error. Also depicted 
is the ground trace to describe the three different satellites 
(Pegasus* Transit and GEOS-B). 

Table 1-A indicates the orbital elements for the simulated satel- 
lites. The error free data used for this simulation were generated for 
each ship location shown. Bias and random errors were added to the data 
as shown in Table 1 1 -A for each tracker and to each ship. The only 
tracker error that was adjusted was Range Rate (R dot). 

Satellite trajectory errors resulted from geopotential errors and 
tracker errors and appear in Table III-A. The errors are shown in the 
conventional H, C, L coordinate system. 

Note that the errors in the satellite position are about the same 
for Range only data as for Range rate (R dot). However* the orbital 
errors are much larger for Pegasus than for Transit or GEOS-B. It will 
be seen later that the orbital errors in Pegasus are the dominant source 

of error in the survey adjustment. 

Table IV-A shows the latitude and longitude for each ship for the 

simulated data obtained for each of the three satellites. The maximum 
elevation (PCA) is shown for each ship also. When the satellite pass 
was almost overhead with high maximum elevation (PCA) poor survey results 
were obtained for each case. 

Table V-A shows the speed and course (heading) for each ship. When 
the survey adjustment was made for position only, an error of 0.5 knots 
and (-0.5) degree heading was always present in the a priori offset for 
the ship. 

128 ' 



Table Vla-A - A geodetic survey adjustment was attempted for each ship 


simulation as shown in Table Vla-A. In one case no solution could be 
obtained. Results are very poor for ship No. 2 which has a high PCA 
and poor geometry. All of the other results in Table Vla-A are 
essentially the same and ship heading has little influence. The errors 
can be largely explained by large orbital errors in the trajectory of 
Pegasus. 

Table VIb-A shows the results for Transit. Ship Mo. 1 gave poor 
results because of high max elevation (PCA). It appears that range 
rate (R-dot) gives less accurate results than does range-only data. 

Also, the velocity of the ship can be adjusted when range data is used. 

If the ship velocity error is large (0.5 knots) it appears to improve 
the survey adjustment when ship velocity is also adjusted. The residual 
error is shown for the ship but the residual does not appear to represent 
the magnitude of the survey error. Also, the R-dot data are sensitive 
to the heading of the ship during a pass. A large effort should be 
made to ascertain a better description of survey error associated with 
range rate (R-dot) and ship heading. The survey errors obtained here 
for R-dot indicate a survey bias could very easily occur. The errors 
are not considered random for the R-dot data. 

Table VIc-A shows the simulated survey adjustment for GEOS-B. Results 
for ship No. 1 are poor because of high elevation (PCA). The R-dot 
data are inferior to the range only data. Better results are also 
obtained for range-only data when ship's speed is adjusted. 
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Table VIIa-A shows the results for range only data to compare results 


with and without orbital errors. (Refer back to Table Vla-A). As 
already indicated 8 Pegasus results are not very good here and require 
further study. Results for ships 1, 3 and 4 are good for the case of 
no orbital errors. However, it may prove difficult to obtain a 
perfect orbit in actual practice. 

Table VIIb-A shows Range-only solutions for Transit. Very little 
difference is indicated between passes with and without orbital errors. 

It should be recalled that the orbital errors were small. In every 
case except high elevation passes (ship No. 1), the results are very 
good. Survey errors are typically 200 feet when ship velocity errors 
of 0.5 knots are present. There is also some indication that heading 
does influence the results but further study is needed to ascertain 
how much. 

Table VIIc-A indicates the simulated results for GEOS-B with and 
without orbital errors. The orbital errors in this case do not strongly 
influence the survey error when ship's velocity is not adjusted. It 
appears feasible to adjust for ship's velocity but in actual practice 
it should be recalled that horizon-to-horizon tracking data are required. 
More effort is required to determine the best procedure to use. However, 
the author would prefer to rely on independent means to determine the 
ship's motion. The ship's ASPS ranging data may prove to be the best 
source for ship's motion. This has not been completely resolved at 
AFETR. 
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TABLE I - A 


ORBITAL ELEMENTS FOR SIMULATED SATELLITES 



PEGASUS 

TRANSIT 

GEOS-B 

Perigee (rum.) 

256 

602 

602 

. . Apogee (rum. ) 

352 

605 

861 

Inclination (deg. ) 

31.7 

90.0 

105.8 


'TABLE II - A 
TRACKER ERRORS 



SITE # 1 

SITE #2 


SHIPS 

R - Bias (ft. ) 

30 

30 

30 

30 

1 R - Bias (ft. /sec.)* 

100 

100 

100 

100 

1 R - Random (ft.) 

25 

25 

25 

10 

1 R - Random (ft. /sec.) 

.15 

.15 

.15 

.10 


* This represents the error remaining pr io r to adjustment. The 

error remaining in R after adjustment was of the order of .1 ft/sec. 
No adjustment was made for R-bias. 
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TABLE III - A 

SATELLITE TRAJECTORY ERRORS IN VICINITY OF SHIPS 



TABLE IV - A 


MAXIMUM ELEVATION OF EACH SATELLITE 
FROM EACH SHIP POSITION 
































TABLE V " A 


TRUE SHIP SPEED AND HEADING* 



*A priori off-sets from those values wore (-.5) knots in speed and 
(-. 5 ) decree in heading for ail simulation runs. 


** Heading is the angle from the reference direction to the direction 
of the ship velocity vector measured positively clockwise as seen 
'by a~ observer looking downward on the ocean surface. 
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TABLE VI b - A 

SHIP POSITION ERROR AT PCA IN SITUATIONS INVOLVING ORBITAL ERRORS 
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TA3LE VI I b - A 

SHIP POSITION ERROR AT PCA IN SITUATIONS WITH AND WITHOUT ORBITAL ERRORS 

(Only Radar Range Data Used in Solutions) 
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TABLE VIII - A 


’ MAXIMUM ELEVATION OF EACH SATELLITE 
FROM EACH LAND BASED SITE 



SITE ■ 

. 

SITE POSITION 

MAX. ELEVATION ] 

SATELLITE 

NUMBER 

N. LAT. (DEG.) 

E. LONG. (DEG.) 

(DEGREES) | 

Pegasus 

1 

28.226 

279.401 

8 . 8 


2 

26.636 

281.732 

13.2 


3 

21.473 

288.868 

42.4 

j 

Transit 

1 

28.226 

279.401 

3 9.5 j 


2 

26.636 

281.732 

4 6.5 1 

© 

3 

21.473 

288.868 

7 8.8 

GEOS-B 

1 

28.226 

279.401 

35, S 


2 

26.636 

281.732 

4 3.1 


3 

21.473 

288.868 

7 6.3 | 






CONCLUSION 


The geodetic location of a ship at sea can be accurately 
determined from pulse radar range-only data. This study has not 
included the relative location of one transponder to another as 
this was considered to be a different problem which involves the 
ASPS measurements obtained by the ship. Also, the ASPS measure- 
ments can degrade the attained accuracy of the transponder survey, 
if not properly used. 

The actual survey accuracy achieved and the number of separate 
satellite passes required depends to a large extent upon (1) the 
accuracy of the reference orbit, (2) the accuracy of the ship's 
tracking data and ship's motion determination and (3) the geometry 
of the different passes used. Where unknown and unmodeled errors 
exist, it will be essential to obtain sufficient passes to examine 
the errors present for each test. It is always desirable to obtain 
satellite passes on each side of the ship for which the geodetic 
survey is desired. 
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ABSTRACT 


"Geodetic Survey Improvement with C-Band Radar” 


AUTHOR 

D. R. Pfingsten - RCA International Service 

Corporation, MTP 


Survey adjustments and improved accuracies were obtained for 
four AFETR downrange Mistram/MRS sites. This presentation 
describes the role played by C-Band tracking radars in obtain- 
ing these results. Comparisons between the C-Band radar 
results and other systems are also discussed. 
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GEODETIC SURVEY IMPROVEMENT WITH C-BAND RADAR 


By; D. R. Pfingsten 


SLIDE 1 


This presentation! describes the role played by C-band track- 
ing radars in reducing survey uncertainties for the Mistram/ 
MRS System. 

This special effort was undertaken because survey uncer- 
tainties at the Mis tram/ MR 3 sites at Grand Turk, Antigua, 
Bermuda, and Trinidad affected the Minuteman III (MM III) 
velocity accuracies after burnout. The flight interval of 
interest was between 300-90 seconds. This accuracy re- 
quirement caused the AFETR to deploy the Mistr am/MRS vans 
at new locations.. It had been shown that if a survey im- 
provement could be effected for the above four X-band sites 
relative to the Florida mainland, this would be a signifi- 
cant contribution towards improving the MM III velocity 
accuracies. 

Briefly, as shown in the figure, the Mistram/MRS system is 
a CW X-band system which consists of transmitter and receiver 
systems at Valkaria and Eleuthera, the rate stations at Miami 
Beach, Pier Road,, Bermuda, Grand Turk, Antigua, and Trinidad. 
The Valkaria and .Eleuthera systems provide unambiguous target 
position both for real time and postflight purposes. 

These rate stations provide a measurement of range sum from 
the transmitter to the missile to the receiver, but the mea- 
surements are ambiguous up to the initialization constant, 
and are used in postflight reductions. 

The Mistram/MRS data are processed postflight in the NITE 
(N Interval Trajectory Estimation Program) which is a minimum 
variance trajectory and error model coefficient estimation 
program designed to process tracking data from free fall and 
powered flight trajectories. All the solutions discussed, 
processed at AFETR, were processed through the NITE Program. 
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SLIDE 2 


In addition to the desirability of the reduction in survey 
uncertainties, data from the first MM III Test 1721 which 
occurred on 16 August 1968 indicated some type of problem in 
the data. One aspect of the problem is shown by the residuals 
from two of the rate stations — Antigua and Trinidad. This 
figure will be shown again. At this time, note that the mag- 
nitude of the residuals, especially late in flight, were 
relatively large. Past experience yielded residual standard 
errors typically on the order of 0.04 ft/sec. 


SLIDE 3 


Theoretical studies performed early in the calibration sat- 
ellite effort had indicated that, even with a moderate amount 
of data, a significant improvement over existing uncertainties 
could be obtained. The limitation on studies of this type is 
that the physical situation assumed in the study may not com- 
pletely describe the actual situation under which the data 
are collected relative to data availability, noise in the data, 
and the significant parameters in the solution. The last two 
limitations apply to the error propagation obtained from the 
actual reduction. 

This slide shows the satellite revolutions for one of the 
two solutions from the GEOS-B data. This particular solution 
involved 9 revolutions. This solution included those revolu- 
tions for which ballistic camera data were obtained. When 
processed with radar only data, this solution is identified 
as NITE 1R. When processed v/ith both radar and camera data, 
it is identified as NITE 6. 

The pulse radars principally used were 0.18, 19.18, 7.18, 

91.18 and 67.18. The first two are located in the Cape 
Kennedy area at Patrick Air Force Base and Merritt Island. 

The last three are located at Grand Turk, Antigua and Bermuda. 

A particular revolution typically included data from three 
to four radars. The ballistic cameras included 1000 mm cam- 
eras at the above sites as well as Homestead Florida, 

Atlantic Field North Carolina, and Trinidad. Data from 
600 mm cameras at Vero Beach Florida and Grand Bahama Island 
were also used as well as some additional data supplied by 
NASA. . ■ 



SLIDE 4 


This slide shows the satellite revolutions for the second 
solution from the GEOS-B data. This particular solution 
involved only radar data. In addition to the GEOS-B data 
shown on the last two slides, passive ballistic camera data 
from the Echo II and Pageos Satellites were collected during 
the same time period the GEOS-B ballistic camera data were 
collected. These data with other available data from past 
programs were used to generate another solution. Finally, 
Mistram/MRS data from three Minuteman III tests were used. 

The failure to significantly change geometry between tests 
and the attendant flame, aspect angle, and staging problems 
associated with powered flight are limitations for these 
data. However, these tests provided data with vehicle motion 
in an east-west direction contrasted to the north-south sat- 
ellite data as well as data from another type of instrumen- 
tation. 


SLIDE 5 

When the decision was made in September 1968 to concentrate 
on the Mistram/MRS Survey Improvement Study, it was also 
decided to confine the radar data used in the solutions to 
range measurements. This decision was based upon studies 
which indicated that exclusion of the radar angle data had 
little effect on the results. It also greatly reduced the 
quantity of data used in the reduction. 

The corrections made to the data prior to input to NITE 
wer e : 

1. Zero-set (including pulse width) 

2. Editing 

3. Transit time 

4. Refraction 

5. Signal strength 

6. Filtering and reduced sample rate 

The zero-set correction was based upon the pre- and post- 
calibrations and included a nominal correction for the 
difference between the calibrate and beacon pulse widths. 
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SLIDE 5 (Continued) 


The refraction correction was computed from a ground index 
value of ref ractivity , determined from measurements, and an 
assumed exponential variation with height. A ray trace com- 
putation was performed. Data below 10° elevation were not 
used to minimize residual refraction errors. The signal 
strength correction was a small correction which attempted 
to account for the variation of beacon delay about its nom- 
inal value. The filtering took the data at 10 samples per 
second through a 101 point filter constrained to pass a 3rd 
degree polynomial. Output data were used at 1 sample every 
10 seconds. A moderate amount of noise reduction was 
achieved. 


SLIDE 6 


The GEOS-B NITE solutions were structured as follows: 

Survey adjustments were made for Grand Turk, Bermuda, and 
Antigua whenever pulse radar data were available. Survey 
adjustments were also made for Trinidad in NITE 6 using 
the available ballistic camera data. Short arc orbital con- 
straints were used in all solutions and the orbital param- 
eters were estimated for each revolution as well as a distinct 
range zero-set for each radar. For all range measurements, 

-7 "" 

a scale factor equal to 5 x 10 was propagated as an un- 
modeled error and was also combined with the a priori noise 
estimates to weight the measurements. A 10° elevation angle 
cutoff was used on all radar measurements to minimize the 
residual refraction effect. In all revolutions involving 
the Bermuda radar, a timing offset, distinct for each revolu- 
tion, was estimated. The original basis for using this 
rather extensive error model was to reduce the risk of obtain- 
ing biased survey adjustments at the expense of less precise 
estimates. For example the inclusion of the Bermuda timing off 
set was not oa» positive evidence that such an offset existed 
but rather on the possibility that it might exist, since it 
alone of all the radars did not receive ETR timing. 
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SLIDE 7 


This slide shows the RMS radar range residual values from 
the two solutions. As can be seen the residuals were small. 
When considered across a single revolution the RMS residual 
values varied from 2 to 6 feet. The residuals were also un- 
biased and the mean value was consistently under 1 foot. 

This of course was not entirely unexpected to the extent 
that range zero-sets were modeled. However, the residuals 
at least gave no positive indication of unaccounted system- 
atic error. 


SL IDE 8 

This slide shows the adjustments obtained. The results are 
presented in terms of a local coordinate system at the bal- 
listic camera site with X-east, Y-north, and Z-up. These 
coordinates are used as a means of presentation. Within the 
RITE Program geodetic coordinates are employed. For the 
islands the radar positions were input relative to the bal- 
listic camera sites and this relative position was held fixed. 
Recall that solutions NITE 4 and NITE 1R are radar only 
solutions, and that NITE 6 is NITE 1R supplemented with bal- 
listic camera data. The pooled results are based upon NITE 4, 
NITE 6 and also upon additional camera and Mistram/MRS data. 
There was no C-band pulse radar at Trinidad. The Trinidad 
adjustment under NITE 1R is a combination of the Antigua ad- 
justment with- an adjustment of Trinidad relative to Antigua. 

As can be seen, the radar solutions agreed very well at all 
sites in the X and Z coordinate adjustments. In these 
coordinates both radar solutions also show good agreement 
with the pooled results. The disagreement in the Y coordinates 
is attributed at least in part to a confounding of these co- 
ordinate adjustments with the Bermuda timing adjustments. 

This effect was also evident in terms of little variance re- 
duction in the Bermuda X and Y adjustments. 


SLIDE 9 

This slide shows the a priori estimates of site uncertainties 
and the uncertainties based upon the pooled results. The 
pooled results are also repeated. As can be seen relative to 
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SLIDE 9 (Continued) 

the a priori uncertainties, the major shifts occurred at 
Trinidad, Antigua X, and Grand Turk Y. The relatively 
large Trinidad shift was explained by an accidental error 
that had occurred in obtaining the previously held values. 

The pooled results were obtained by computing a weighted 
mean across the various solutions. In this computation 
which involved pooling results from various types of instru- 
mentation and geometry, the tacit assumption was made that 
the relative weights within a solution, based on the computer 
outputs, were correct but scaled improperly. 


SLIDES 10 & 11 


These slides show a comparison between the MRS residuals 
shown earlier and the residuals from the same test with the 
new survey values. It can be seen that the residuals are 
reduced particularly late in flight. The&e residuals are 
of course only one aspect of the problem. For this particu- 
lar test the survey shifts produced a significant shift in 
the trajectory estimated by the system. The shifted trajec- 
tory proved to be more consistent internally and also agreed 
better with trajectory solutions from other systems. The 
large shift at Trinidad accounts for a major part of these 
changes. While this large shift was indicated from other 
sources including the Mistram/MRS data, the similar results 
from the GEGS-B data were important. This was due to the 
fact that the test shown was the first test supported at the 
Trinidad site by the MRS van and it would have been difficult 
to establish by other means that the shift determined was in 
fact a survey error and was not some other unknown error 
functionally similar to survey. 


SLIDE 12 


Subsequent to the adoption of the values previously shown, 
four additional Minuteman III tests provided Mistram/MRS 
data suitable for adjustments. In these solutions an accurate 
trajectory was the primary objective and the survey adjust- 
ments were performed as a means to properly weight the 
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SLIPS 12 (Continued) 


observations for trajectory determinations. This slide 
shows the mean values from these four tests and the un- 
certainties as of August 1969. These mean values corres- 
pond to an additional shift relative to the shifts previously 
shown on Slide 9. The fact that they are in general small 
and within the uncertainties quoted supports the consistency 
of the changes adopted. The one exception to this is the 
Bermuda X shift. The problem here does not, however, appear 
to be due to the GEOS-B data used. 


SLIDE 13 


This study was the first attempt at AFETR to use C-band 
pulse radar data from satellites in short arc solutions for 
survey improvement. The conclusions reached were as follows: 

1. The method can properly determine any major inconsis- 
tencies. 


2. Weaknesses in the adjustment will be identifiable from 
internal comparisons. 

3. When the results were internally consistent, compari- 
son with adjustments obtained by other means were for 
the most part consistent within the uncertainties quoted. 


4. The results were obtained relatively fast compared to 
results using ballistic cameras. 
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TEUCTURE OF NITE SOLUTIONS 
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ADM RESIDUAL 
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ACOUSTIC DATA REDUCTION AND CALIBRATION FOR THE 
GEOS -II C-BAND RADAR MARINE GEODESY EXPERIMENT 
IN THE BAHAMAS 


by 

A. G. Mourad, A. T. Hopper, D. M. Fubara, and G. T* Ruck 

Battelle Memorial Institute 
Columbus Laboratories 


INTRODUCTION 

In late 1969, Battelle was involved in a NASA experiment in the Bahamas 
the overall objective of which was to determine the feasibility of using C-band 
radar for ship positioning and for establishing marine geodetic control points. 

It was the first such experiment initiated by NASA and was aimed ultimately at 
the practical application of existing satellite and space technology for better 
understanding of the earth, particularly the complex ocean environment. Battelle ' s 
initial role was to handle that portion of the experiment requiring acoustic 
techniques . 

Results obtained, using the Apollo ship Vanguard and the GEOS-II 
satellite, were excellent, despite existance of some gaps in the data obtained and 
the fact that general-purpose acoustic system components were used. For example, 
standard errors of about 0.24 arc second and three meters were obtained in determining 
the latitudes and longitudes, and depth of the three ocean-bottom transponders 
respectively. 

Figure 1 shows the location of the experiment and indicates, generally, 
how the ship position was determined relative to ocean-bottom-mounted acoustic 
transponders during observed passes of the GEOS-II satellite. The new acoustic 
procedures which were developed, referred to here as Transponder Location by 
Surface Positioning (TLSP), permit use of any type of accurate surface- ship- 
positioning data available; they involve measurement of acoustic slant ranges 
from the ship to any number of underwater acoustic transponders or hydrophones 
and consideration of the velocity of sound. TLSP can utilize all available 
acoustic data and lends itself to a weighted least-squares determination of 
transponder coordinates. 
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FIGURE 1. LOCATION OF C-BAND RADAR MARINE GEODESY EXPERIMENT 


3 


The feasibility of establishing marine geodetic control points on the 

ocean floor using bottom-mounted acoustic transponders had previously been 

demonstrated in November, 1968, by the Battelle/Industry Group experiment in 

( 6 ) 

the Pacific Ocean . A standard point error of ± 15 to ± 20 m was obtained in 
determining the coordinates of the control point that was established about 120 
miles west of Los Angeles using a LQRAC airborne DME line-crossing technique. 
Although the airborne DME system is capable of providing fairly high geodetic 
accuracy at sea, it is range- limited (up to a few hundred miles). Therefore, 
the use of satellite systems is attractive because of the world-wide coverage 
possible by their use. 

Present indications are that if high-accuracy in marine geodetic measure 
ment can be achieved, many difficult problems will be well on the way to solution 
Such problems include, for example, determining ocean circulation, establishing 
reference points for ocean tidal measurements, determining ocean- floor spreading, 
providing "ground truth" for calibration of satellite altimetry and evaluating 
shipboard positioning systems, making accurate gravity measurements, determining 
geoidal heights, and making deflection-of-the-vertical measurements. 

The use of C-band radar for marine measurements is desirable for two 
reasons: (1) C-band radar can be operated in a ranging mode (measure distances 
from ship to satellites) which provides a particular advantage in that errors 
in ship velocity are not as critical as they are with the Doppler satellite 
techniques (the only satellite system operational at sea) and (2) several ships 
already have C-band radar systems on board. 

GENERAL DESCRIPTION OF THE EXPERIMENT 


Figure 1 shows the general location of the experiment conducted in 
the Bahamas. This location was chosen: (1) because of its proximity to the test 
area for the Vanguard (the experiment was a task added to the primary Vanguard 
mission) and (2) because of the availability of the four underwater acoustic 
transponders previously established by the U. S. Navy Compass Island ship in 
about 16,000 feet of water. The four transponders formed a quadrangle having 
about 3.5 mile sides. However, during the experiment it was possible to collect 
data from only three of them. 
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The experiment involved six major elements: (1) the Apollo ship 

’’Vanguard 5 ’ and its instrumentation systems, (2) the GEOS-II satellite, 

(3) three underwater acoustic transponders, (4) a C~band radar land-based 
tracking network, (5) a LORAC surface positioning network, and (6) the Navy 
Transit satellites. The ship instrumentation systems used include: (1) the 

FPS-16 C-band radar, (2) the SRN-9 Doppler satellite navigation receiver, 

(3) the Inertial/Star Tracker System (INS), a LORAC receiver, (4) the Bathymetric 
Navigation System (BNS) , and (5) timing and computer systems. 

Basically, simultaneous measurements were made of (1) C-band radar 
ranges (between the ship and the GEOS-II satellite), (2) acoustic ranges (between 
ship and three underwater acoustic transponders, (3) surface LORAC positions 
(between ship and three land LORAC transmitters, and (4) the Inertial/Star Tracker 

coordinates . 

The ship was positioned continuously over the triad formed by the three 
ocean-bottom-mounted acoustic transponders during the GEOS-II and Doppler satellite 
passes. Acoustic data were collected during 5 of 7 GEOS-II passes, 22 of 32 
Doppler satellite passes, and 11 acoustic calibration runs. 

DESCRIPTION OF ACOUSTIC SYSTEM USED 

The acoustic system used consisted of two major components: (1) the 

three ocean-bottom-mounted acoustic transponders and (2) the shipboard Bathymetric 

Navigation System (BNS). The BNS system on board the Vanguard is capable of 

interrogating ocean-bottom-mounted transponders, receiving their replies, and 

recording the acoustic slant ranges to the transponders. This capability of the 

(1 2 ) 

BNS was utilized in this experiment. Operation of the BNS is as follows ’ 

The BNS transmits an interrogation signal with a frequency of 16.0 KHz ± 50 Hz 
and a pulse width of 15 ± 1 milliseconds (ms) to the transponders. The leading 
edge of the interrogation signal activates a digital counter which starts counting 
the pulses from an 800- pps clock. The reply signal from a transponder is applied 
to a signal processing unit which recognizes the center of the returned pulse 
and stops the counter. The time interval measured by the counter represents the 
range from the BNS transducer to the acoustic transponder (in yards, an 
acoustic velocity of 4800 ft/sec is assumed) plus several fixed delays. To 
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correct the time interval and obtain the true range, the BNS subtracts 175 ms 
(8 ms for transponder- turn-around time, 160 ms for BNS signal-processing time, 
and 7 ms for signal-recognition time) from the measured time interval. 


ACOUSTIC TECHNIQUES FOR DETERMINING 
GEOMETRY OF TRANSPONDERS 


Two techniques have been developed for determining the coordinates 
and orientation of the ocean-bottom transponders. These are: 

(1) Acoustic line-crossing technique 

(2) Ocean-surface positioning technique. 

The basic line-crossing technique involves the measurements of: 

(1) the ship's heading, (2) the acoustic slant ranges or the two-way travel 
times of an acoustic signal from a ship to two adjoining ocean-bottom trans- 
ponders, and (3) the velocity of sound^^ . The ship, allowed to travel at a 
constant speed and heading, crosses the lines joining the transponders and records 
the points of closest approach (PCA) to each of the transponders and the acoustic 
slant ranges versus time during the crossings. The sum of two slant ranges 
is minimum when the ship is in the vertical plane containing the two transponders. 

The minimum sums can be obtained by least-squares curve fitting or graphically 
by plotting the horizontal ranges against time. From the ship's heading and 
course, and minimum distances, forward and backward azimuth, the relative coordinates 
of the transponders are determined. The technique is sensitive to ship heading 
and must depend on good data, particularly during the PCA's. Batte lie's experience 
in the Pacific, as well as in the Bahamas, indicates that adverse operational 
conditions and/or equipment malfunction can occur, thus limiting the value of 
using such techniques. 

The ocean-surface positioning technique, known as transponder location 
by surface positioning (TLSP), involves the use of a computer code incorporating 
any available accurate surface-ship coordinates (e.g., LORAC, SINS, DECCA) , 
acoustic slant ranges from the ship to the transponders, and velocity-of- sound 
profile to determine the positions of the transponders. Basically, each transponder 
was independently positioned relative to a set of LORAC surface positions. The 
LORAC positions were given in terms of geodetic latitude and longitude. These 
positions were transformed to a space rectangular coordinate system X, Y, Z where 
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the Z-axis coincides with the polar axis of the reference ellipsoid and the 
X- and Y-axes are situated in the equatorial plane, perpendicular to each other, 
and X-axis going through the O-meridian. The original Navy coordinates for the 
transponders were used as approximate coordinates (X , Y q , z q )> and an observation 
equation V = AX + L, was written for each selected LORAC position. The computed 
slant ranges SR q are given by: 

SR = [(x - x.) 2 + (y - y ) 2 + (z - z.) 2 ] 1 ^ 2 

o o 1 o i o i 


where 


(x, z^) = coordinates of ship converted from LORAC geodetic 

positions at times i; i = 1, n where n is the 

total number of observations. 


The adjustment was carried out with the following matrices: 


X = approximate values of unknowns (in this case they 
o 

represent the approximate transponder coordinates) 

X = adjusted values of unknowns 
a 


X 

SR b 


L 



= observed slant ranges as determined from the two-way 
acoustic travel times and velocity of sound (all slant 
ranges were corrected for refraction): 



dSR 
o 

S (Xo, Y o> Z o> 


The least-squares solution of the observation equations, 
V = AX + L 

is given by 

X = -N _1 U 

whe re 


172 



1 


N = A ! PA 
U = A'PL 

P weight matrix. 

After the adjusted X, Y, Z coordinates of the transponders were 
obtained, their orientation was determined by computing the inverse geodetic 
problem. The solution gives also the variance covariance matrix and the 
eigenvectors and eigenvalues of this matrix which are used in the error ellipsoid 
computations . 


ACOUSTIC TECHNIQUES FOR DETERMINING SHIP POSITIONS 

Ship positions relative to ocean-bottom-mounted acoustic transponders 

( 6 ) 

can be determined by either of two techniques developed at Battelle' and for 
which computer codes exist. 

The first technique which employs the least-squares solution for a 
two-dimensional computer model is based on the intersection of three circles 
in a plane. The Z-coordinate is held fixed at Z = 0. One transponder is used 
as the origin of the local coordinate system with the positive X-axis correspond- 
ing to east and the positive Y-axis corresponding to north. 

An observation equation AX + L = V is written and solved in matrix 
form. The general solution is similar to that already described. The unknowns 
in this case are the rectangular ship coordinates (X, Y) at times of the acoustic 
readings . 

The second technique, which employs the three-dimensional computer 
code is based on the intersection of three slant ranges corrected for velocity 
of sound and refraction. The coordinate system is an earth-centered XYZ system 
as in TLSP. The general equation for the slant ranges from the ship to the 
transponders is given by: 

(X - X.) 2 + (Y - Y.) 2 + (Z - Z.) = SR 2 . 
ill 


The solution of the quadratic equation results in two positions. These two 
positions will result in two values for X. The desired position is that 
corresponding to the larger X value. 
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Veloc ity-of-Sound Measurement 

Actual measurement or determination of velocity of sound during the 
fine of operation is essential for accurate reduction of acoustic data. At the 
time of the experiment, the Vanguard did not have a capability for making velocity- 
of-sound measurements. Therefore, it was necessary that existing data be 

(2 

utilized. Two ve locity-of-sound curves were provided for the area of operation ’ 
Figure 2 shows the curve made by the Compass Island velocimeter during July, 1969, 
with a mean value of 4989 ft/sec and the other curve represents the mean of 
five observations taken in the area during the months of August and September 
in different years prior to 1967 with a mean value of 4964 ft/sec . These 
mean-velocity values were not used in the final computations. Instead, data 
points were read directly from the two velocity-of-sound curves into the ray- 
tracing program. 


DISCUSSION OF RESULTS 


The results of TLSP are summarized in Tables 1, 2, 3, 4, 5, and 6. 

Tables 1 and 2 give the adjusted geodetic coordinates of the transponders and their 
standard errors in both curvilinear coordinates (based on the I960 Fischer 
Ellipsoid) and cartesian coordinates. The data used to obtain these coordinates 
included approximately 1500 acoustic ranges from the ship to the three underwater 
acoustic transponders taken during the acoustic calibration runs and the GEOS-II 
satellite passes. The ship positions consisted of LORAC surface positions and 
SINS positions during these runs. Several corrections were made to the raw 
data in an attempt to eliminate systematic errors. The corrections were for 
ship velocity effects on the acoustic data, antenna offset from the BNS transducer, 
and corrections for velocity of sound and acoustic refraction. 

The least-squares solution of the observational data included weighting 
criteria. Table 3 shows the six different cases that were applied for determining 
the proper weighting criterion for the experiment data. 

Table 4 shows the coordinates of the transponders and their standard 
errors as determined by using the two available sound velocity profiles. It 
also gives the transponder coordinates as were determined by the Navy. 
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FIGURE 2. SOUND VELOCITY CURVES FOR BAHAMAS EXPERIMENT 
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TABLE 1. ADJUSTED GEODETIC COORDINATES (CURVILINEAR) BASED 

ON NAVY SOUND VELOCITY PROFILE (1960 FISCHER ELLIPSOID) 


Transponder 
No , 

0 

Latitude - 
(North) 

i 

' 9 

ill 

Longitude 

(West) 

0 t 

- X 

II 

Height 

(m) 

Standard Errors 

V"> CT x (,,) CT h (m) 

2 

27 

8 

14.75 

76 

23 

14.60 

-5018.0 

0.24 


3.38 

3 

27 

8 

22.05 

76 

20 

03.73 

-5016.5 

0.23 


2.65 

4 

27 

5 

44.92 

76 

21 

35.51 

-5016.1 

0.19 

0.20 

2.99 



TABLE 2. ADJUSTED 
ON NAVY 

GEODETIC COORDINATES (CARTESIAN) BASED 
SOUND VELOCITY PROFILE 



Transponder 

No. 

X 

Y 

Z 

Standard 

Errors 

(m) 

(m) 

(m) 

cr (m) 

A 

CTy(m) 

a (m) 
z 

2 

1,335,766.98 

-5,516,086.75 

2,889,502.18 

4.74 

3.32 

7.76 

3 

1,340,846.95 

-5,514,749.83 

2,889,702.74 

5.46 

1.96 

7.07 

4 

1,338,912.43 

-5,517,486.64 

2,885,401.47 

5.56 

5.07 

4.37 







TABLE 3- CRITERIA FOR WEIGHTING OBSERVATIONS 
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TABLE 4. TRANSPONDER COORDINATES AS DETERMINED 

BY USING DIFFERENT SOUND VELOCITY PROFILES 
(1960 FISCHER ELLIPSOID) 


TLSP Determinations Using 

Transponder Navy Navy Sound Std. Average Sound Std. 

Numbe r Determination Velocity Profile Error Velocity Profile Error 


2 

9 

27 

o 

00 

16.32" 

27 

00 

o 

14.75" 

0.24" 

27 

00 

o 

14.43" 

0.24" 


X 

76 

23' 

16.56" 

76 

23' 

14.60" 

0.19" 

76 

23' 

13.93" 

0.19" 


h* 

- 

4966 

.7 m 

- 

5018.0 m 

3.38 m 


- 4998 

• 3 m 

3.38 m 

3 

9 

27 

08' 

21.54" 

27 

08' 

22.05" 

0.23" 

27 

08' 

21.21" 

0.24" 


X 

76 

20' 

03.42" 

76 

20' 

03.73" 

0.19" 

76 

20' 

04.13" 

0.21" 


h* 

- 

5023 

.4 m 

- 

5016 

.6 m 

2.76 m 


-4998. 

7 m 

2.89 m 

4- 

CD 

27 

05' 

44.16" 

27 

05' 

44.92" 

0.19" 

27 

05' 

45.45" 

0.19" 


X 

76 

21' 

38.28" 

76 

21' 

35.51" 

0.20" 

76 

21' 

35.44" 

0.20" 


h* 

- 

4977 

.6 m 

- 

5016 

.6 

2.99 m 


-4995. 

2 m 

2.98 m 


* A geoidal undulation of 60 m (i.e., geoid is 60 m below the ellipsoid) was assumed. 
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Although the standard errors for each transponder did not change much, influence 
of sound velocity variations on the coordinates was more noticeable,, particularly 
in the depth determination. This influence on the depth was of the order of 18 
to 21 m for the three transponders. Furthermore, the transponder ellipsoidal 
depths as determined by TLSP are about the same for all three transponders while 
they were different in the Navy determination. 

Table 5 shows clearly these differences in the three transponder depths 
(ocean surface to bottom) as determined by TLSP and the Navy. The largest 
difference between TLSP, using Navy velocity, and the Navy determination is 
about 51 m in Transponder 2. Transponder 3 has the smallest difference of 7 m. 

Table 6 presents the results of the error ellipsoid computations for 
all three transponders. The largest errors are associated with the depth values 
for Transponder 2 and 3. Both these transponders did not have sufficient ship- 
track data over them. Transponder 4 has one ship track data over it, which is 
indicated perhaps in the small error component in the depth axis. In general 
the geometry of the error ellipsoids indicate that better results could be 
obtained if data all around the transponders were taken rather than the one- 
sided data used here. 

Preliminary reduction of the acoustic data taken during the calibration 
runs revealed many obviously erroneous data at the "critical points", such as at 
the times of PCA's and line-crossings. Therefore, only four "best” runs (Numbers 
8, 9, 10, and 3A) were selected and reduced. Figure 3 shows a plot of the 
acoustic calibration runs as determined by surface LORAC positions. The results 
are summarized in Table 7 and Table 8. Because of the inaccuracy of the data in 
line crossing determinations, it was necessary to develop the TLSP technique. 

Table 9 presents primarily the results of the distance and azimuth 
determinations by the TLSP and line crossing techniques. In addition, the 
distances between the transponders and their forward and backward azimuths were 
also computed from the Navy given coordinates. 

Table 10 summarizes the results of the angles between the three 
transponders from TLSP, line crossing techniques, and the Navy computed values. 

Once the final coordinates were determined for the three ocean bottom 
transponders, the relative ship positions during the GEOS passes are then deter- 
mined by using program TRISPR. The results of the ship positions during GEOS 
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FIGURE 3. LORAC SHIP TRACKS DURING ACOUSTIC CALIBRATION RUNS 
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TABLE S. TRANSPONDER DEPTHS 


Transponder 

Number 

Navy 

Determination 

Bates lie Determinations from TLSF 
Using Sound Velocity Profiles by 

Navy Average. 

2 

4906.7 m 

4958.0 

4938.3 

3 

4963-4 m 

4956,6 

4938,7 

4 

4917-6 m 

4956.6 

4935-2 



TABLE 6. PARAMETERS OF CONFIDENCE 
ADJUSTED COORDINATES AT 

: ELLIPSOID FOR 
95% CONFIDENCE 

THE 

LEVEL* 

( 8 ) 


Transponder 

Number 

Semi-axis 
in Meters 

Axes Orientation Angles with Respect 
to the Positive Directions of the 
Geodetic Coordinate Axes 
X Y 2 

2 

15.30 

37° 

12' 

123° 

56* 

103° 

28"* 


3.51 

54° 

27’ 

35° 

38' 

87° 

42' 


22.44 

80° 

23' 

99° 

37 ' 

13° 

40' 

3 

15.05 

15° 

46* 

85° 

24' 

105° 

02 ' 


3.70 

91° 

32 ' 

11 ° 

47 * 

78° 

19* 


20.72 

74° 

18’ 

100 ° 

49' 

IS 0 

13® 

4 

15.13 

21 ° 

32' 

116° 

16' 

82° 

16* 


18.47 

64° 

30* 

46° 

55' 

126° 

10 s 


3.34 

80° 

14 s 

54° 

28' 

52 c 

44 s 


* This is not the conventional standard error ellipsoid (whose confidence 
level is 20 %) but a confidence ellipsoid scaled such that the probability 
that the point being determined lies within it is 95%,. (8) 
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TABLE 7. TRANSPONDER DISTANCES DETERMINED FROM LINE CROSSINGS 


Run No. 

T 2 “ T 3 
<m) 

t 2 “ t 4 

(m) 

T 3 - ' I 4 
(m) 

8 

5182.2 

5498.9* 

- 

9 

- 

5488.8* 

5393.8* 

10 

5428.5 

5325.0 

5428.5 


(5457.4) 

(5405.6) 


3A 

5323.0 

- 

5498.6 


(5320.6) 


(5519.3) 

Avg. 

5320.8 

5437.6 

5440.3 


(5315.7) 

(5451) 

(5443.7) 


* Poor data, not graphed 

( ) Distance determined from individual values at 

times of crossings 


TABLE 8. TRANSPONDER AZIMUTHS DETERMINED FROM LINE CROSSINGS 


Run No . 

t 2 - t 3 

(o) 

T 3 - t 2 
(») 

T2 - T 4 
(°) 

\:> T2 

t 3 - t 4 
( °) 

W 3 

(°) 

8 

89 

269 

146 

326 

- 

- 

9 

- 

- 

152 

332 

207 

27 

10 

- 

- 

151 

331 

207 

27 

3A 

85 

265 

- 

- 

201 

21 

Avg. 

87 

n 

149.6 

329.6 

205 

25 
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TABLE 9 . COMPARISON OF TLSP WITH LINE-CROSSINGS AND NAVY RESULTS 


Dist., in 
Meters 

Az., degrees Lines TLSP Line Crossings Navy 


Dist. 

V T 3 


5260 

.9 m 


5315 

.7 m 


5321 

.0 tn 

Dist. 

V T 4 


5358 

. 6 m 


5451 

tn 


5409 

.2 tn 

Dist . 

V T 4 


5457 

. 1 tn 


5443 

. 7 tn 


5503 

.5 m 

Az. 

T 2 -T 3 

87° 

32' 

17.06" 

87° 

00’ 

00" 

88° 

15' 

27=07" 

Az. 

V T 2 

267° 

33’ 

44.12" 

267° 

00’ 

00" 

268° 

16' 

55.17" 

Az . 

V T 4 

149° 

22' 

29.55" 

149° 

36' 

00" 

149° 

57* 

56.30" 

Az . 

T 4- T 2 

329° 

23' 

14.72" 

329° 

36’ 

00" 

329° 

58" 

41.09" 

Az . 

V T 4 

207° 

36' 

02.26" 

205° 

00' 

00" 

208° 

20' 

51.76" 

Az . 

T 4- T 3 

27° 

35’ 

20.43" 

25° 

00' 

00" 

28° 

20’ 

08.52" 


TABLE 10. COMPARISON OF TRANSPONDER ANGLES 


Angle 


TLSP 



Line Cross 

ing 


Navy 


o 

S 

5! 

O 

f 

tl 

0 

t 

iT" 

a 2 

61 

50 

12.49 

62 

36 


61 

42 

29.23 

ff 3 

59 

57 

41.86 

62 

00 

B 

59 

56 

03.41 

ff 4 

58 

12 

05.71 

54 

24 

00 

58 

21 

27.43 
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passes are shown graphically in Figure 4 at one-minute interval. Figure 5 shows 
the same tracks as determined by the LORAC surface ship positions. 

It is significant to note that the ellipsoidal height, h, of the 
ship (initially assumed to be-60 meters) was recovered accurately in the computed 
acoustic time history of the ship motion. The maximum variation in the computed 
heights was about 9 m. The mean value was 60.05 m with a standard deviation of 
± 1.65 m. These results indicate the confidence in the adjusted coordinates of 
the transponders and their derived depths. The fact that the computed ship 
positions were derived from only three intersecting ranges (no redundancy) but 
corrected ray paths indicate the reliability of both TLSP and the acoustic ray 
tracing programs. 


SUMMARY OF RESULTS AND CONCLUSIONS 


An improved technique for determining the geometry and orientation of 
underwater acoustic transponders has been completed, and analysis of data obtained 
during the Bahamas experiment indicates that the technique will make possible 
more accurate measurements at sea than other techniques previously known available 

The results of TLSP adjusted geodetic coordinates in both curvilinear 
and cartesian coordinates are shown in Tables 1 and 2. Analysis of TLSP results 
revealed that standard errors of less than 0.24 arc seconds were obtained in the 
determination of the latitudes and longitudes for all three transponders. 

Moreover, the standard error of depth determination was about 3 meters at 
depths of about 5000 meters. Such remarkable accuracy, particularly in depth 
determination will make it possible to establish vertical reference datum and 
provide, for example, "ground truth" for future satellite altimetry measurements. 

About 1500 observation equations were used with nine unknowns (trans- 
ponder coordinates) in a least-squares solution to determine the transponder 
coordinates and orientation. Corrections applied to the data include velocity 
of sound and refraction corrections, antenna offset correction between the LORAC 
antenna and the BNS transducer and correction for ship speed. In addition, from a 
thorough investigation of several possible weighting criteria coupled with 
statistical tests, proper weights were applied to the observation in the least- 
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FIGURE 4. ACOUSTIC SHIP TRACKS DURING GEOS PASSES 


27° 09' N 


20 



DURING GEOS PASSES 




■ 21 


(8 9) 

squares solution' ’ , Two velocity-of-sound curves were used. The results 

indicate sensitivity of the TLSP technique to the velocity of sound curves used . 
Velocity of sound effects were most noticeable in depth determination. Results 
proved to be sensitive to ship speed, and antenna offset correction made a negligible 
effect on the precision estimates of the adjusted coordinates. As expected,, 
the increased number of observations improved the precision of the adjusted 
coordinates of the transponders. 

The line-crossing technique was applied to the data first but preliminary 
data reduction revealed there were many obviously erroneous data recorded at the 
points of closest approach (PCA’s) and at the times of line crossing. Although 
the line-crossing technique is relatively simple with respect to data reduction, 
it is operationally time consuming. Moreover, it is apparent that the technique 
is sensitive to ship heading and course, and its effectiveness depends on the 
obtaining of reliable and accurate data, particularly during PCA's. 

The time-history of ship positions with respect to the adjusted 
transponder coordinates as determined by TLSP proved the accuracy and reliability 
of the TLSP and the ray tracing programs. The ellipsoidal heights of the ship 
positions with respect to the transponders were recovered with standard deviation 
of ± 1 . 65 m. 

On the basis of the results of this experiment and results of an earlier, 
similar experiment, several conclusions can be drawn and a number of recommendations 
can be made. The results achieved to data show how the newly developed techniques 
for data reduction and analysis can improve the accuracy and potential use of data, 
even when they are not as precise as desired. Furthermore, all past investigations 
have involved the use of off-the-shelf equipment designed for general-purpose 
use and not specifically for geodetic programs. Indications are that with 
appropriate limited modification of present equipment, the potential for achieving 
high accuracy at sea is substantial and that eventually investigating difficult 
problems such as ocean spreading will be quite possible. Accordingly, the following 
steps are recommended: 

(1) Develop an operational test area suitable for conducting controlled- 
condition experiments for testing hardware and newly developed, 
techniques 

(2) Modify and improve present hardware designs, particularly these for 
underwater acoustic transponders and shipboard receiving equipment. 
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I. INTRODUCTION 


One of the continuing projects at AFETR is the maintaining and 
understanding of instrumentation accuracy. This is accomplished 
both by hardware modifications and by software analysis. The spe- 
cific aspect that this report is concerned with is the calibration 
Information that can be obtained for the C-band radar by collecting 
and analyzing radar data from satellites with C-band beacons (soft- 
ware analysis). 

AFETR's interest in instrumentation calibration utilizing satel- 
lites was first proposed in 1957. A satellite calibration project was 
approved in 1960 under AFETR Project 8938. When all formal satel- 
lite projects were transferred to ESD by Headquarters AFSC in 1962, 
the new calibration task was ESD Project 5930, Task 5930. 03. Under 
this project, an AFETR calibration satellite was formulated [ l] and 
was unsuccessful in achieving orbit on 13 July 1966. Instrumentation 
calibration utilizing satellites was continued by use of the Agena 10 
and 11 as well as the NASA GEOS-B (See Reference [ 2] ). This par- 
ticular study was accomplished utilizing data obtained from GEOS Bf 

There are many aspects to radar calibration. Idealistically, we 
would like to obtain a permanent set of calibration coefficients that 
would always be applied to the radar data and whose values would 
change very slowly. In this way, the permanent set of calibration 
values would only have to be updated periodically. However, a more 
realistic situation is that some of the important calibration values are 
only constant over short periods of time and, therefore, need to be 
re-evaluated on a timely basis. For example, the azimuth (A) and 
elevation (E) zero set biases may remain constant over a few weeks 
or in some instances may only remain constant for one or two days. 

In a multiple radar solution, the A and E zero set biases would not 
be important since range (R) basically determines the solution and A 
and E can b e easily self-calibrated. However, for a single radar 

* See references [ 4 ] and [5] . 
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solution, A and E are very important and cannot be easily self-calibrated. 
Therefore, the calibration of A and E zero set biases constitute a signifi- 
cant contribution to improving the accuracy of single radar solutions. 
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II. ANALYSIS 


In order to obtain a good calibration of A and E biases it is necessary 
to exercise the radar over a range of A and E values. That is, a dynamic 
calibration is certainly more valuable than a static calibration. The tar- 
get board static calibrations are useful for obtaining the gross A and E 
biases, but a dynamic calibration is needed for the fine grain correction 
needed to obtain the ultimate accuracy. Therefore, a satellite with a 
C-band beacon can be invaluable for obtaining good dynamic angle cali- 
brations for the radar system. 

ETR is in the unique position of being able to immediately process 
radar data on the CDC-3600 at Cape Kennedy due to the real time data 
transmission link between the MIPIR radars and the computer. It is, 
therefore, possible to immediately obtain timely angle calibration values 
from a multiple radar satellite calibration pass to be applied to radars 
when they are tracking a different vehicle. The CDC-3600 computer pro- 
gram, which is capable of producing an orbital fit for many radars, is 
the GPTP program (See reference [ 10] ). If the same processing is done 
post flight where the radar data is carefully edited and corrected, then 
the NITE program [ 3 ] (run on the IBM 7094) is utilized. A comparison 
of the main features between GPTP and NITE is shown in Table 1. 

The purpose of this report is to evaluate the usefulness of the GPTP 
program in producing good near real time calibration values. The cri- 
teria is based on a comparison of these values with the post flight reduction 
done on the NITE program, which is assumed to be the most accurate data 
processing technique for obtaining angle calibration values. The set of 
data used for these comparisons were obtained utilizing GEOS B in 
March 1968 and May 1968. A detailed description of the ETR Radar 
Calibration using GEOS B is given in reference [ 11] . 

For the March 1968 experiment, radar data were collected and ana- 
lyzed over five GEOS B revs of data: 801, 807, 808, 813, 814. The radar 
coverages is shown in Table 2 and specified details of data processing can 
be found in references [ 6 ] and [ 7 ] . The comparison between GPTP and 



TABLE 1 

COMPARISON OF GPTP AND NITE SOLUTIONS 
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TABLE 3 

GPTP VS NITE FOR MARCH 1968 DATA 
(Angles in Milliradians) 

















NITE for A and E angle biases only are given in Table 3 for radars 0. 18, 

19. 18 and 7. 18. The omitted comparisons were either due to obvious 
hardwartff problems or data not collected for that particular revolution. 

The specific details of this March experiment were documented in 
a Secret report [ 7 ] . The results presented in this report constitute the 
unclassified part of the report. From Table 3, it is seen that there is an 
excellent agreement between the GPTP and NITE solutions for angle biases 
where the great majority of the azimuth bias differences are < . 04 mr. As 
noted in [ 7 ] , there was some hardware discussion concerning reasons why 
some of the zero set levels changed. For example, for radar 0. 18 there 
was a malfunctioning parametric amplifier over revs 801-808 and there 
was a possibility of a mislevel error. With the exception of the 0. 18 bias 
levels, which changed significantly probably due to a hardware problem, 
the other bias levels for 19. 18 and 7. 18 appear fairly stable, at least in 
the sign of the correction. There seems to be a fair amount of consis- 
tency between revs 801 and 807 and between revs 808 to 814 for radar 
19. 18. It appears that angle bias results may be stable for one-half to 
one day. 

Another GEOS B experiment was run in May 1968 (see references [ 8 j , 

[ 9] ) where we also obtained comparison computer runs on GPTP and NITE. 
In this particular comparison, the zero set biases were assumed to be con- 
stant from revs 1672 through 1700, or over a period of two days. The 
radar coverage for these revs is given in Table 2. Both the GPTP and 
NITE programs were run under the constraint of constant biases. Table 4 
shows the results of the May 1968 comparison. 


TABLE 4 - GPTP vs NITE FOR MAY 1968 DATA 
Includes Revs 1672, 1673, 1674, 1685, 1686, 1687, 1698, 1699 and 1700 . 

(Angles in Milliradians) 
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As in the March 1968 experiment (Table 1), there is good agreement 
between the GPTP and NITE angle zero set biases,, The differences 
between the NITE and GPTP angle biases are mostly <. 04 mr with the 
exception of 91. 18 elevation difference of . 09 mr and 12. 18 azimuth 
difference of . 06 mr . Although these differences are larger than those 
found in Table 1, they are still small enough in order to utilize the GPTP 
results as meaningful. The main reason for the additional variation is 
probably due to the fact that the biases were assumed to be constant 
over two days. It seems to be apparent that more consistent results 
can be obtained if we restrict our bias evaluation to less than two days. 

An analogy would be if we performed the standard radar pre-calibration 
two days before the radar was to be used. Obviously, for those tests 
that would require good angle calibrations, the dynamic satellite cali- 
bration should be done as close as possible to the time the radar will 
be used for observing test data. More explicit information concerning 
the May 1968 experiment can be found in Reference [ 9 ] . 

Summing up, it can be stated from the data comparisons shown in 
Tables 3 and 4 that the angle calibrations from the GPTP program are 
usable in a near real time sense to obtain improved radar test data. It 
is perfectly conceivable that we could obtain updated angle biases from 
a multiple radar GPTP solution and apply the angle bias calibrations to 
a single radar test - all in the same day. 
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ABSTRACT 


C-band radar tracking of the GEOS-B satellite 
by the Vanguard Apollo tracking ship has been used, 
in conjunction with land-based C-band radar tracking, 
to obtain estimates of the ship position. For a two 
revolution dockside solution, a position agreeing with 
the surveyed position to within 17 meters was obtained. 
At sea radar tracks by the Vanguard were reduced using 
the ship inertial navigation system (SINS) data for 
the change in ship position after the beginning of 
track. Both one and two revolution solutions were 
used for ship positioning. In all cases, poor fits of 
the radar data to the orbit determined by the land- 
based stations gave indications of problems with the 
SINS data and station recoveries in error by 2 km or 
greater. Comparison with a LORAC position for one case 
also showed discrepancies on the order of 3 km. 



SECTION 1.0 
INTRODUCTION 


The precise determination of land-based geodetic 
positions through the use of satellite tracking is a 
well established art. In general, the procedure is to 
use minimum variance statistical estimation techniques 
such that a set (or sets) of orbital parameters and sets 
of tracking station coordinates provide a "best fit" of 
the tracking measurements to the satellite orbit. The 
coordinates of some stations must, of course, be con- 
strained in some manner (such as fixing latitude and 
longitude) in order to obtain a non-singular solution. 

A similar procedure may also be used to estimate 
the geodetic location of a tracking instrument located 
on a ship. It is, however, necessary to have information 
about the motion of the ship during the tracking period. 
This information may take several forms, including: 

(1) inertial navigation data, (2) ranges from the ship 
to fixed acoustic transponders, and (3) radio fixes to 
land-based transmitters as via the LORAC network. The 
results presented in this paper use primarily the first 
type of information. For some of the tests discussed, 
some data of both of the other types does exist and 
will be used in data reductions yet to be completed. 

With perfect knowledge of the change in ship 
position from one measurement point to the next, the 
accuracy of station position estimation can be estab- 
lished through the use of an error analysis simulation, 
or through the use of ship data when the ship is known 
not to be in motion. The latter technique was actually 
used. This has the advantage of gaining familiarity 
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with the use of actual data and an assessment of the 
errors which may exist in the ship-borne tracker. 

With the ship at sea, the technique used was to 
calculate the change in ship position from its position 
at a reference time (generally at the beginning of a 
satellite track) using the ship inertial navigation 
system data. The Wallops A/Omega program, which accepts 
data from a moving tracker, was then used to estimate 
the position of the ship at the reference time in a 
minimum variance solution which also determined the 
orbital parameters. The measurement weighting used, 
however, was such that the orbit was determined primarily 
by land-based tracking stations. 



SECTION 2.0 
SELECTION OF DATA 


Although C-band radar data almost always includes 
azimuth and elevation measurements, such data is of much 
lower accuracy than are range measurements . In an 
orbit determination using data from more than one radar, 
it is possible to utilize range measurements only, 
with no degradation of overall orbit accuracy. For 
this reason, the orbits for ship positioning were ob- 
tained using land-based radar range measurements. In 
most cases, data from only two land-based locations was 
available. However, error analyses have shown that 
this is not a serious limitation so long as the land- 
based stations are not widely separated from the ship 
which is being positioned. 

The FPS-16 type radar on the Vanguard makes 
angle measurements just as its land-based counterparts 
do. However, since calibration and stabilization 
problems are greatly more difficult on the ship than 
on land, the angle data would be expected to be, and is, 
far less accurate in terms of positioning accuracy than 
is the range data. The estimation of a two (or three) 
dimensional ship position using range data alone from 
a single satellite pass can be done, in most cases, but 
is rather sensitive to errors of almost any nature. 

Two satellite passes provide a much greater positioning 
capability as should be expected simply on the basis 
of the extended exercise of geometry. (For GEOS-B, 
e.g., tracking of the satellite on successive revolutions 
automatically means a pass east of the ship and a pass 
west of the ship.) For this reason, considerable em- 
phasis has been placed on obtaining ship positions using 
two passes of satellite track. 
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The only data which has thus far been available 
for obtaining ship movement, relative to some fixed 
position or time, has been from the ship inertial 
navigation system. In the data reductions, it was 
assumed that the SINS position might contain an error 
at the reference time, but that the error remained 
constant throughout the one or two satellite passes. 

It was recognized in advance that this was only an 
approximation, and perhaps a poor one. The analysis 
reported in this paper had as one of its objectives 
to determine the degree of validity of this approxima- 
tion. 


It should be noted that the data reductions herein 
reported were made using "raw" C-band range data from the 
ship, with the reduction performed in a moving tracker 
program. Computer programs on board the ship are de- 
signed to provide C-band data referenced to a fixed 
point. This data, generally called low-speed data, 
contains, in addition to the SINS errors, the effects 
of angle biases and any radar range biases which may 
exist. Since this transformation introduces large 
errors from the biases alone, it was decided to use 
the radar data in its raw form so that radar errors 
themselves, if they exist, could be more clearly 
identified. 

One such error source known to exist in the C-band 
data is a range bias. Far stronger position recoveries 
would be possible if the radar data were sufficiently well 
calibrated for its bias to be negligible. Unfortunately, 
calibration procedures used to date have not satisfied 
this requirement. For all data reductions performed to 
date, it has been found necessary to require the data re- 
duction process to estimate an independent radar bias on 
each pass. 


206 



SECTION 3.0 
DOCKSIDE TESTS 

Using two GEOS-B tracks by the Vanguard C-band 
radar with the ship in Port Canaveral, it has been possible 
to check out the radar and data handling and reduction 
processes under near ideal conditions. Lack of motion by 
the ship removed the complications of the moving tracker. 

And the position of the tracker was tied to a well surveyed 
point. Position recovery under these conditions should 
give an indication of the accuracy that should be expected 
when the ship motion is perfectly accounted for. 

Table 1 shows the radar position recovered when tracks 
of Revolutions 7972 and 7973 were used in conjunction w r ith 
land-based tracking by Wallops and Bermuda radars. The 
estimated position, when compared to the survey position also 
shown in Table 1, shows a difference of only 17 meters. This 
figure should be considered a measure of the accuracy of the 
estimated position relative to Wallops and Bermuda. 

For both radar passes, independent radar range biases 
were estimated along with the ship position. Biases of -67 
and -17 meters were recovered, indicating both significant 
biases and biases which differ from pass to pass. 

The ship radar and land based radar residuals (differ- 
ences between the observed ranges and the calculated ranges 
to the fitted orbit) for Revolutions 7972 and 7973 are shown 
in Figures 1 and 2, respectively. The raw radar range residuals 
are denoted by ISTA19. The residuals denoted by ISTA66 refer 
to the 1 orv speed data transformed by the ship computer. The 
trends in these residuals are an indication that something 
was wrong in the transformation process. The ISTA19 residuals 
seem quite comparable to the land-based radar range residuals. 
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TABLE 1 


DOCKSIDE SHIP TEST 
GEOS-B REVS 7972 - 7973 
NAD-27 POSITIONS 


Adjusted Survey- 


Latitude 

28° 

24 ’ 

31.8" 

28° 24’ 

31.4 

Longitude 

279° 

23 ' 

44.1" 

279° 23 ' 

44.5 

Height 


7 . 2 

m 

14.6 

m 


Orbit Determination by: 

Wallops FPS-16 
Bermuda FPQ-6 
Bermuda FPS-16 


20 8 




igure 1, Residuals for GEOS-B Rev 797 
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Figure 2. Residuals for GEOS-B Rev 7973 



SECTION 4.0 
AT SEA TESTS 


At sea Vanguard tracks on three different days have been 
reduced, in 'all cases using the SINS data for relating the 
position of the ship during track to its position at the 
beginning of track. When tracks on successive revolutions were 
available, the initial position used was that at the beginning 
of the first track. 

Results for the three different tests will first be 
given. Some analysis will then be made of the recovered 
ship positions and the implications of the large measurement 
residuals obtained for all at sea tests. 

It should be noted that a height recovery was not 
attempted for any of the at sea tests. This constraint was 
used because of the limited amount of data on some passes 
and the necessity for recovery of range biases on ail passes. 

In addition, the probability of height improvement did not 
seem very high. 

4.1 ONE PASS SOLUTION 

Table 2 gives the estimated ship position at the 
beginning of track on GEOS-B Revolution 7991. The only- 
available comparison position is the SINS position, also 
listed in Table 2. Differences are approximately 440 meters 
in latitude and 30 meters in longitude. A reliable 
estimate for the accuracy of the SINS position is, however, 
not available. 

The residuals for the ship range measurements after 
the ship position (and radar bias) estimation are shown in 
Figure 3. Residuals for the land based radars were negligibly 
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TABLE 2 


AT SEA SHIP TEST 
GEOS-B REV 7991 


Latitude Longitude 


SINS 

28° 23’ 

Ad jus ted 

28° 23' 


Orbit Determination by: 

Wallops 

Wallops 

Bermuda 


22.4" 

280° 15’ 

Y — 1 

00 

o 

36.7" 

280° 15’ 

16.8 


FPQ-6 

FPS-16 

FPS-16 



Ship Range Residuals for GEOS-B REV 7991 



(sroroui) s^rnpisotf oSuny 


dn 


small on the scale of this figure and were consequently not 
plotted. The explanation for the large systematic residuals 
almost certainly lies in errors in the SINS relative posi- 
tioning inaccuracy. A preliminary analysis of the residuals 
assuming the simplest error models for a well damped SINS 
system, indicates that such models can do only a very poor 
job of explaining the actual residuals. A simple velocity 
error alone must be in excess of 3m/sec to give the best 
residual fit. This would indicate rather serious SINS 
problems during this pass of an unknown nature. 

4.2 TWO PASS SOLUTIONS 

Two sets of two pass solutions have been reduced, 
with the results shown in Tables 3 and 4. For one of these 
tests, a LORAC position is available for comparison with 
the estimated position and the SINS position. If the LORAC 
position is accepted. Table 3 says that the adjusted posi- 
tion is in error by about 500 m in longitude and 3000 m 
in latitude. These differences appear to be largely 
attributable to SINS errors, as discussed below. 

Figure 4 shows the ship range residuals after the 
ship position (and range bias) estimation on Revs. 8003-8004 
Figure 5 shows the . corresponding residuals for Revs. 8010- 
8011 . For both tests, land based radar range residuals are 
essentially negligible on the scale of the graph and for 
that reason are not plotted. 

As might be expected, SINS systematic errors over a 
period of two hours do not necessarily take a very simple 
form and it is not possible to account easily for the residu 
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TABLE 3 


AT SEA SHIP TEST 
GEOS-B REVS 8003 - 8004 


Latitude Longitude 


SINS 

27° 

7* 

17.5" 

283° 

39 ' 

0.6 

Ad jus ted 

27° 

9 ' 

12.9" 

283° 

39 * 

10.9 

LORAC 

27° 

7 ' 

32.0" 

283° 

38 ' 

53.0 


Orbit Determination by: 

Bermuda FPS-16 
Antigua FPQ-6 


TABLE 4 


AT SEA SHIP TEST 
GEOS-B REVS 8010 - 8011 


Latitude Longitude 


SINS 

27° 

6 ' 

41.8" 

283° 

38 ' 

10.5 

Adj us ted 

27° 

7' 

26.5" 

283° 

37' 

38.0 


Orbit Determination by: 

Wallops FPS-16 


Bermuda FPS-16 



shown in either Figures 4 or 5. For Figure 4, however, a 
fair fit to the residuals is obtained by assuming velocity 
errors of 0.3 m/sec S and 0.1 m/sec W. These errors corres- 
pond to position errors of about -400 m in longitude and 
about 1600 nj in latitude. These numbers correspond to 
about half the observed LORAC -A/Omega differences for 
latitude - and in the right direction - but are in the 
wrong direction for improving longitude. 

It may be noted that the SINS position shows less 
agreement in Table 3 with the estimated position than it 
did in Table 2, probably due mainly to the longer time 
span. However, the SINS velocity errors necessary to 
account for observed residuals are lower by about an order 
of magnitude. 

Results for Revs. 8010- 8011 appear to be somevvhat 
worse than those for Revs. 8003-8004. The SINS and esti- 
mated positions of Table 4 show closer agreement, but the 
residuals shown in Figure 5 are larger. Most of the data 
taken fro this test is on Revolution 8011. Reduction of 
this pass by itself (and land based data) shows residuals 
of the same magnitude and shape as was found for Rev. 7991, 
giving again indications of rather poor SINS performance. 
Accordingly, further analysis of this test was postponed 
until additional data is available. When available, acoustic 
transponder data will be combined with the radar data to 
obtain an improved set of ship positions. 


4.3 AT SEA DATA ANALYSIS 

The results of much of the analysis of the sea test 
data have been given above. In general, the large resi- 
duals obtained for the ship radar are simply not explicable 
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3S iduals for GEOS-B REVS 8003 8 800. 
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on any basis other than errors in the SINS data used for 
the ship motion history. Several of the simpler SINS error 
model forms have been modeled in a simulation program to 
see their effects on the recovered ship position and on 
the measurement residuals. With a very limited amount of 
data,, as was the case for all tests, it is impossible to 
solve for a large number of error model coefficients. 

This difficulty is due to the fact that a number of the 
error model terms have similar forms over the periods for 
which data is available. In such a situation, one or two 
parameters may be chosen to effectively absorb the error 
in all parameters. In all the analysis made to date, no 
more than two components of velocity error have been use- 
ful in accounting for observed residuals. 

The large errors in the estimated ship positions 
can, as has been emphasized, be attributed to large SINS 
errors, with no accounting made in the data reduction for 
such errors. Some improvement can be made, in most cases, 
by the solution for one or two selected parameters. As 
has been noted, however, the magnitude of SINS errors 
appears to differ greatly from day to day. This suggests 
that the data might be potentially much better than it 
generally is. 

The velocity errors which have been quoted as 
accounting for observed residuals should not be interpreted 
as being actual velocity errors, particularly for the two 
revolution solutions. The velocity error best accounting 
for the 8003-8004 residuals was stated to be 0 . 33/ 19 7 ° 

The speed and heading of the ship was completely changed 
between Revolutions 8003 and 8004, with some expected 
perturbations on the SINS system. The velocity errors 
would then be expected to be different on the two passes. 
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SECTION 5.0 
CONCLUSIONS 


The results presented in the preceding sections 
have shown that C-band radar data from the Vanguard Is 
potentially useful for accurate geodetic positioning 
provided the relative ship movement can be accurately 
accounted for. The results obtained and the implications 
for further analysis may be summarized as follows: 

1. A two satellite pass position estimation 
accuracy on the order of 15 meters is possible 
using a ship C-band radar when relative ship 
position between data points is known from 
other sources. 

2. Ship position estimation using SINS data for 
relative ship motion gives position errors 
which may exceed 3 km. (One single pass 
result, GEOS-B revolution 8004, gave answers 
which residual analysis indicated should be 
accurate to 500 meters.) 

3. Apparent SINS errors observed are much greater 
than quoted SINS errors. 

4. Modeling of SINS errors over a one pass period 
appears possible, with the solution for speed 
and heading errors. 

5. The crucial item in ship positioning is an inde- 
pendent source of relative position. The use 

of acoustic transponders for this purpose should 
be one perfectly adequate solution. 
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